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1 . Introduction 

The  last  five  years  have  brought  considerable  advances  in 
the  theory  of  the  boundary-val ue  problem  of  physical  geodesy  in  the 
formulation  of  Molodensky,  which  is  the  determination  of  the 
physical  earth's  surface  from  gravity.  These  advances  have  been 
accomplished  mainly  through  the  work  of  T.  Krarup,  L.  Hbrmander 
and  F.  Sansd. 

The  present  report  is  devoted  to  a review  of  this  work.  Its 
aim  is  to  introduce  the  reader  to  the  basic  ideas  and  geodetically 
important  results,  which  are  sometimes  hidden  between  formidable 
mathematical  technicalities.  We  shall  thus  attempt  what  mathema- 
ticians call  a "heuristic  exposition",  for  mathematical  details 
the  reader  will  be  referred  to  the  original  papers.  The  treat- 
ment of  the  linear  problem  in  gravity  space  in  sec. 8 is  new. 

The  problem  of  Molodensky  may  be  formulated  briefly  as 
follows;  given,  at  all  points  of  the  physical  earth's  surface  S , 
the  gravity  potential  W and  the  gravity  vector  £ , to  deter- 

mine the  surface  S . The  potential  W can  be  determined  by 
leveling  combined  with  gravity  measurements;  this  gives  the 
potential  up  to  an  unknown  constant  which,  however,  can  be  found 
indirectly  by  other  methods,  especially  distance  measurements. 

The  magnitude  of  the  gravity  vector  £ , which  is  gravity  g , is 
measured  by  gravimetry,  and  the  direction  of  £ , which  is  the 
plumb  line,  is  obtained  by  astronomical  measurements  of  latitude 
4>  and  longitude  A . It  is  assumed  that  these  measurements  have 
been  corrected  for  luni -solar  tidal  effects  and  other  temporal 
variations,  so  that  our  problem  is  independent  of  time.  We  further 
suppose  that  the  effect  of  the  atmosphere  has  been  taken  into 
account  by  appropriate  reduction.  Hence,  the  space  outside  the 
surface  S can  be  considered  as  empty. 
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We  thus  assume  that  the  earth  is  a rigid  body  which  rotates 
with  constant  and  known  angular  velocity  oj  around  a fixed  axis, 
which  passes  through  the  earth's  enter  of  mass.  This  center  of 
mass  will  be  taken  as  the  origin  0 of  a cartesian  coordinate 
system,  the  x^  axis  coinciding  with  the  axis  of  rotation. 

The  gravitational  potential  V is  a harmonic  function 
outside  S . For  large  values  of  the  radius  vector 


r = xi 


2 2 2 
Xi  + X2  + X3 


(1-1) 


it  has  an  expansion  in  spherical  harmonics  of  the  form 


GM  Y(0,A) 

Hi)  = ^ + —2 + + ...  . 


(1-2) 


where  G is  the  gravitational  constant,  M denotes  the  total 
mass  of  the  earth,  and  Y^(0,x)  are  Laplace  surface  harmonics, 

0 (polar  distance)  and  x (longitude)  forming  together  with  the 
radius  vector  r a system  of  spherical  coordinates  related  to 
the  cartesian  coordinates  x = (x^,  x^.  x^)  by 

x^  = r sin0  cosx  , 

X2  = r sin0  sinx  , (1-3) 

x^  = r COS0 

The  condition  that  the  coordinate  origin  0 coincides 
with  the  center  of  mass  implies  that  the  spherical  harmonics  of 
first  degree  vanish  identically: 
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Yj(e.A)  E 0 . (1-4) 

SO  that  V must  have  the  form 

V(x)  = — + O(^)  for  r - - . (1-5) 

The  gravity  potential  W is  then  given  by 

W(x)  = V(x)  + + x^)  . (1-6) 

It  will  also  be  assumed  that  the  surface  S is  a one-to- 
one  image  of  the  sphere  and  that  it  is  a smooth  surface,  being 
differentiable  as  often  as  required. 

It  may  be  questioned  whether  Molodensky's  problem  thus 
formulated  is  to-day  geodetically  relevant  at  all.  On  the  one 
hand,  the  prerequisites  for  Molodensky's  problem,  especially 
continuous  coverage  of  the  whole  earth's  surface  by  gravity 
measurements,  are  still  far  from  being  realized;  on  the  other 
hand,  there  are  many  more  date  of  different  kind,  such  as  satellite 
data,  that  transcend  the  frame  of  Molodensky's  problem  and  must  be 
handled  by  data  combination  techniques  such  as  least-squares 
col  1 ocati on . 

To  these  questions  we  may  answer  as  follows.  From  a 
practical  point  of  view,  the  integral  formulas  arising  in  the 
solution  of  boundary- val ue  problems  are  often  computationally 
more  convenient  than  collocation  and  retain  their  importance  if 
gravity  data  are  available  to  a sufficient  extent,  at  least  locally 
(cf.  Mor i tz , 1975) . From  a theoretical  point  of  view,  the  geodetic 
boundary-value  problem  represents  an  especially  interesting  and 
significant  special  case,  whose  importance  for  the  conceptual 
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structure  of  geodesy,  from  the  time  of  Clairaut  to  the  present  day 
can  hardly  be  overestimated  (curiously  enough,  the  theory  was 
always  far  ahead  of  the  data  available  at  the  time).  In  fact,  the 
consecutive  stages  in  the  development  of  the  boundary- val ue  problem  - 
Clairaut,  Stokes,  Molodensky  --  always  served  as  measures  of 
perfection  for  geodetic  theory  and  set  new  standards. 

Even  today  Molodensky's  problem  is  not  yet  completely 
clarified  from  a mathematical  point  of  view,  with  respect  to 
existence  and  uniqueness  of  the  solution,  in  spite  of  the  decisive 
progress  made  in  the  last  few  years;  it  remains  a challenge  to 
theoreti ci ans . 

Let  us  now  try  to  get  a first  grasp  of  the  mathematical 
nature  of  Molodensky's  problem. 

The  gravity  vector  £ can  be  expressed  in  terms  of  measured 
gravity  g and  of  astronomical  latitude  $ and  longitude  A as 


r g COSO  cosA  1 

£ = ' g COSO  s 1 nA 
g s i no 


(1-7) 


In  space  the  vector  £ and  the  potential  W may  be  considered 
functions  of  the  rectangular  coordinates: 

1 = ><2’  ^^3)  ’ W " W(x^,  x^,  x^)  (1-8) 

On  the  earth's  surface  S , they  are  functions  of  two  surface 
coordinates,  for  which  we  may  take  the  astronomical  coordinates 
0 and  A : 


S = i(o,A)  , 


T?  = W(o,A)  ; 


(1-9) 
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the  overbar  denotes  restriction  of  spatial  functions  to  the  surface 
S , whereas  underlining  characterizes  vectors  (and  later,  matrices). 

iNow  £ may  be  expressed,  in  a certain  sense,  as  a function 
of  S and  y , symbolically 

£ = F(S,W)  . (1-10) 

This  means  that,  given  the  surface  S and  the  gravity  potential 
y on  it,  the  gravity  vector  £ on  S is  then  uniquely  determined 
and  can  be  computed. 

In  fact,  this  may  be  done  as  follows.  Let  S and  W be 
given.  Compute  the  centrifugal  potential  on  S (which  can  be  done 
since  the  surface  S is  supposed  to  be  given  and  consequently  the 
coordinates  x^,  x^,  x^  of  the  surface  points  are  known)  and 
subtract  it  from  y ; this  gives  the  gravitational  potential  Y 
on  S . From  Y on  S we  get  the  potential  V outside  S by 
solving  Dirichlet's  boundary  value  problem,  which  has  a unique 
solution.  Now 

£ = grad  V + centrifugal  force 

(grad  denoting  the  gradient)  can  be  computed  outside  S and,  by 
the  continuity  of  first  derivatives,  also  on  S , giving  £ 

Thus  £ is,  in  fact,  uniquely  determined  by  S and  W , so 
that  (1-10)  holds. 

Suppose  now  that  it  were  possible  to  solve  (1-10)  for  S : 

S = $(y,£)  . (1-11) 

This  would  express  the  earth's  surface  S in  ■t''rms  of  W and  £ , 

solving  Molodensky's  problem. 
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This  is  probably  the  conceptually  simplest  formulation  of 
Molodensky's  problem.  However,  the  transition  from  (1-10)  to  (1-11) 
is  mathematically  extremely  difficult.  If  S , TJ  and  ^ were 
simple  real  numbers  and  F were  an  ordinary  function  (supposed 
sufficiently  smooth),  then  the  solution  of  (1-10)  for  S would 
be  straightforward.  The  existence  of  such  a solution  is  guaranteed 
by  the  elementary  implicit  function  theorem. 

In  fact,  however,  the  "function"  F in  (1-10)  is  a rather 
complicated  nonlinear  operator,  and  the  existence  of  a solution 
(1-11)  is  by  no  means  obvious.  There  are  implicit  function  theorems 
for  nonlinear  operators  (e.g.  D i eudonne  , 1960  ; Loomi s and  Sternberg, 
1968  ; Schwa rtz , 1969  ; Sternberg , 1969 ) , but  the  conditions  for  their 
application  are  not  satisfied  in  the  geodetic  case.  It  was  the 
merit  of  Hdrmander  (1975)  to  have  found,  by  a mathematical  tour  de 
force , an  implicit  function  theorem  that  is  applicable  to  the 
geodetic  boundary-val ue  problem. 

To  get  some  first  insight  into  the  matter,  let  us  forget 
all  mathematical  difficulties  and  proceed  formally  as  if  S , W , 
and  £ were  simply  real  numbers  and  F were  a simple  functions. 
Since  ¥ is  given,  it  can  be  considered  fixed  once  and  for  all, 
so  that  (1-10)  becomes  a function  of  S only; 

I = f{S)  . (1-12) 

To  further  simplify  the  notation,  we  write  g instead  of  £ , 

obtaining 

g = f(S)  . (1-13) 

Thus  S is  simply  given  by  the  inverse  function 

S = f-^g) 
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so  that  the  implicit  function  problem  reduces  to  an  inverse  function 
probl em. 

To  practically  find  this  inverse  function,  that  is,  to 
solve  (1-13)  for  S , we  may  apply  the  usual  procedure  for  solving 
nonlinear  equations,  namely  1 inearization. 

Let  us  introduce  an  approximation  to  the  earth's  sur- 

face S and  let  g be  the  corresponding  gravity  vector,  related 
to  by  ( 1- 13 ) : 

g^  = f(SJ  . (1-15) 

Write,  formally. 


S = S + AS 

O 

g = g^  + Ag 


(1-16) 


and  apply  Taylor's  theorem  to  (1-13): 


+ Ag  = f(S^  + AS)  = 

= f (S  ) + f ■ (S  )aS  , 

' o o' 


omitting  quadratic  and  higher  terms.  In  view  of  (1-15)  this  becomes 
Ag  = f ' ('S  )aS  . (1-17) 

o 

The  formal  solution  of  this  equation  is 


AS 


Ag 


(1-18) 


Let  us  link  these  ideas  with  the  conventional  approach  to 
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Molodensky ' s problem.  Here  S is  the  telluroid  and  is  normal 

gravity  on  it;  Ag  is  the  usual  gravity  anomaly  referred  to  the 
earth's  surface  (it  is  here  possible  to  disregard  the  original 
vector  character  of  Ag  and  regard  it  as  a scalar  quantity)  and 
AS  is  represented  by  the  height  anomaly  c characterizing  the 
separation  between  earth's  surface  S and  telluroid  . Thus 

(1-18)  becomes 


MAg 


(1-19) 


where  M = Tf  ' (S  ) I denotes  the  linear  Molodensky  operator 
computing  ? from  Ag  ; practically  one  uses  Stokes'  formula 
with  suitable  corrections. 

Higher  approximations  may  be  obtained  by  Newton's  method. 
Combining  (1-15),  (1-16)  and  (1-18)  we  get 


S 


1 


- )1 

'o' 


(1-20) 


where  we  have  written  S^  instead  of  S to  indicate  that  by  this 
equation  we  get  better  approximation  S^  rather  than  the  true 
value  S itself.  By  repeated  application  of  this  formula  we  get 
successive  better  approximations  S^,  S^ 


^2  = S * 


r - 1 

! ^ • 


S3  = Sz  . (S^)! 


g - f(sj  1 
g - f(s,)] 


(1-21) 


Graphically  Newton's  procedure  is  illustrated  by  Fig.  1. 
The  unknown  abscissa  S for  the  given  ordinate  g is  approached 


Figure  1 


Newton ' s Method 


by  following  the  broken  line  with  arrows. 

The  convergence  of  Newton's  procedure  is  known  to  be  very 
good,  namely  quadratic:  there  is  a constant  K independent  of  n 
such  that 


. - S I < K is  - S , 

' n n- 1 


fl-22 
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The  following  sections  2 to  4 will  deal  with  a detailed 
study  of  the  linearized  problem.  The  remaining  sections  will  be 
devoted  to  two  approaches  to  the  nonlinear  problem,  Hbrmander's 
appraoch  is  logically  straightforward  , using  an  iterative  inverse 
function  technique  basically  similar  to  Newton's  method  but  mathe- 
matically extremely  involvpd.  The  second  approach,  due  to  Sansd, 
first  transforms  the  o*^iginal  free  bounoary  value  problem  (the 
boundary  S is  "free",  that  is,  unknown)  into  a fixed  boundary- 
value  problem  by  means  of  a Legendre  transformation,  thereby 
essentially  reducing  the  mathematical  complexity. 

2 . Krarup's  Linearization 

In  the  usual  linearization  of  Molodensky's  problem,  the 
telluroid  is  introduced  as  the  surface  formed  by  the  set  of  points 
Q such  that  Q lies  on  the  same  ellipsoidal  normal  as  the 
corresponding  point  P at  the  earth's  surface  and  that  the  normal 
potential  U at  Q is  equal  to  the  actual  potential  W at  P ; 
cf.  (Heiskanen  and  Mor i tz , 1967  , p . 292 ) . 

In  his  third  letter  on  Molodensky's  problem  that  was 
circulated  among  the  members  of  the  lAG  Study  Group  on  Mathematical 
Methods  in  Physical  Geodesy  but  unfortunately  never  published, 
Krarup  (1973)  gave  a more  general  formulation  of  the  linearization 
which  is  also  suitable  for  studying  the  nonlinear  problem. 

In  this  mere  general  formulation,  the  tel  1 uroi d ^ is 
now  an  arbitrary  given  surface  close  to  the  earth's  surface  S , 
the  points  Q of  which  are  in  some  one-to-one  correspondence 
with  the  points  P of  S ; cf.Fig.  2.  We  also  introduce  a normal 
potential  U which  constitutes  an  analytic  approximation  to  the 
actual  potential  W ; U is  usually  taken  as  the  gravity 
potential  of  an  equi potential  ellipsoid. 

n 


It  should  be  mentioned  that  the  first  rigorous  formulation  and 
linearization  of  Molodensky's  problem  has  been  given  by  Meissl 
(1971). 


Figure  2 


The  telluroid  I as  an  approximation 
to  the  earth's  surface  S 


Let 


X = grad  U (2-1) 

denote  the  normal  gravity  vector,  in  the  same  way  as 

£ = grad  W (2-2) 

expresses  the  actual  gravity  vector. 

Since  I and  U are  given,  we  can  compute  U and  at 

Q , that  is,  Ug  and  Xq  • potential  W and  gravity  £ are 
supposed  to  6e  given  on  S (in  the  notation  of  Section  1,  they  are 
W and  £ ),  we  know  it  at  every  point  P on  S , that  is,  we 
know  W and  £ . We,  therefore,  can  compute  the  differences 
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called  potential  anomaly  and  (vectorial)  gravity  anomaly,  respectively. 

By  appropriate  definitions  of  the  telluroid  it  is  possible 
to  make  one  of  the  two  quantities  (2-3)  and  (2-4)  equal  to  zero.  In 
the  usual  definition  of  the  telluroid  mentioned  at  the  beginning 
of  this  section,  we  have  = Wp  and  therefore 

AW  = 0 . (2-5) 


In  this  definition,  points  P and  Q are  supposed  to  lie  on  the 
same  ellipsoidal  normal.  Since  the  ellipsoidal  normal  through  P 
is,  strictly  speaking,  not  known,  it  would  be  theoretically  more 
correct  to  define  Q by  the  three  conditions 


U 


(2-6) 


Here  (fi  and  X are  defined  by 

r 1 

; Y COS<()  COSX 

_Y='YCOs<|>sinx  (2-7) 

Y s i ncf) 

L 

in  complete  analogy  to  (1-7);  thus  the  normal  latitude  ()>  and 
longitude  x determine  the  direction  of  the  normal  gravity  vector 
;y  . in  the  same  way  as  and  A define  the  direction  of  £ 

The  surface  formed  by  the  points  Q in  this  manner  has  been  called 
"normal  surface"  in  (Mor i tz , 1964  ) . Krarup  ( 1973)  calls  it  "Marussi 
telluroid"  because  the  three  "Marussi  coordinates"  potential, 
latitude  and  longitude  are  identified. 

In  this  way,  the  potential  anomaly  aW  can  be  made  zero. 
Somewhat  surprising  at  first  sight  is  that  also  the  gravity 
anomaly  A£  can  be  made  to  vanish.  This  requires  defining  the 
points  Q of  the  telluroid  by 
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= £p  • (2-8) 

Expressing  this  vector  condition  in  terms  of  magnitude  and  direction 
of  the  vectors  involved,  we  get  three  conditions 


Y 


Q 


Q 


(2-9) 


which  again  completely  determine  Q . Since  g,  ^ , A may  be  called 
"gravimetric  coordinates",  the  corresponding  locus  of  points  Q 
has  been  called  by  Krarup  the  "gravimetric  telluroid";  for  it,  in 
fact , 

Ai  = 0 . (2-10) 

After  these  possible  specializations,  let  us  return  to  the 
general  case  in  which  both  aW  and  A£  are  nonzero.  As  usual, 
we  define  the  disturbing  potential  T by 

T = W - U , ‘ (2-11) 

W and  U referring  to  the  same  point  (this  distinguishes  T from 
the  potential  anomaly  aW  , in  which  W and  U refer  to  different 
points!). 

On  substituting 

Wp  = Up  + Tp  (2-12) 


we  get  from  (2-3)  and  (2-4) 
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Tp  + Up  - Ug  = AW  , (2-13) 

£p  - Xq  = • (2-14) 

Let  us  now  preceed  with  the  linearization.  We  put 

£ = vector  QP  (2-15) 

(see  Fig.  2)  and  systematically  neglect  all  quantities  of  second 
and  higher  order  in  c . It  is  well  known  and  easy  to  see  that 
quantities  such  as  T and  a^  have  the  same  order  of  magnitude 
as  £ .So  also  T , T^,  etc.  are  quantities  of  second  order  to 
be  neglected. 

By  a Taylor  expansion  restricted  to  linear  terms  we  get 
Up  = U^  + grad  U • c = + y • 1 , (2-16) 

where  the  dot  denotes  the  inner  product  of  two  vectors.  Let  us 
proceed  in  the  same  way  with  the  normal  gravity  vector: 

Xp  = Xq  + grad  X ‘ i • (2-17) 

What  is  grad  x 2 To  see  this,  let  us  write  this  equation  in  index 
notation,  using  the  summation  convention  (summation  over  an  index 
occurri ng  twi ce) : 


P,i 


^ I'll 

8x. 


Y • + M . ^ ^ , 
Qfi  ij  D 


(2-18) 


where 


_ _ 3 /3U  V _ 3^U 

^ij  ^ 9x.  ax.  ^ax. ' ■ ax. ax. 

3 31  ^3 


(2-19) 


Hence  grad  y is  nothing  else  than  the  matrix 


M = M . . i = 
- L 13  J 


ax . ax . 
1 3 


(2-20) 


formed  by  the  second  derivatives  of  the  normal  potential  U ; in 
other  words,  it  is  the  second-order  normal  gradient  tensor.  There- 
fore, we  may  write  (2-17)  as 


Yq  = Yp  - Mg. 


(2-21) 


It  is  clear  that  ^ in  (2-16)  and  M in  (2-21)  refer  to  point 

Q . 

Let  us  similarly  expand  : 


Tp  = + grad  T • ^ . 


Now,  however,  grad  T is  already  small  of  first  order,  so  that 
grad  T • £ is  of  second  order  and,  therefore,  negligible.  Thus, 
consistent  with  our  linear  approximation,  we  simply  have 


• 

P Q 


(2-22) 


The  insertion  of  (2-16),  (2-21),  and  (2-22)  into  (2-13) 
and  (2-14)  now  gives 


Furthermore , 


Ip  ■ Xp  = (gi^aci  W)p  - (grad  U)^ 

= grad  (W  - U)p 
= (grad  T)p 
= (grad  T)^  , 

for  the  same  reason  as  (2-22).  We  thus  finally  get 


T + = AW 

9 

(2-25 

grad  T + M;; 

= , 

(2-26 

in  which  T and  grad  T refer  to  Q 

, as  well 

as 

Y and 

M . W 

have  used  the  matrix 

T 

notation  a b 

for  the  inner 

product 

1 • ^ 

the  transpose  of  a 

being  denoted 

u T 

by  a 

These  two  equations  will  be 

basic  for 

our 

further 

devel op 

ments.  Let  us  solve 

(2-26)  for  £ 

t assuming 

M 

invertible, 

^ = M"^(Aja  - grad  T)  , (2-27 

and  substitute  into  (2-25): 


T + - grad  T)  = aW 
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t 

I 


or 

T - /M'^grad  T = aW  - . (2-28) 

On  putting 

m = - M"\  (2-29) 

we  get 

T + m'^grad  T = aW  + m^A^  (2-30) 

This  equation,  which  holds  on  the  telluroid  I , consti- 
tutes the  fundamental  boundary  condition  for  the  linearized 
Molodensky  problem.  It  is  a generalization  of  the  "fundamental 
equation  of  physical  geodesy"  (Heiskanen  and  Mor i tz , 1967 , p . 86 ) , 
just  as  (2-25)  is  a generalization  of  Bruns'  formula  (ibid.  , p . 8 5 ) . 

Various  Forms  of  the  Boundary  Condition.  - Let  us  intro- 
duce new  coordinates  q^  by 


"l  = S’ 

S - q2(x,.  X2.  X3) 

(2-31) 

S = S'S-  S'  S> 

or  briefly 

(2-32) 

and  let  us  assume  that  the  inverse  transformation 


(2-33) 
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also  exists.  More  specifically,  we  shall  select  to  be  the 

cartesian  components  of  the  normal  gravity  vector: 


(2-34) 


It  is  clear  that  one-to-one  relations  (2-32)  and  (2-33)  exist,  at 
least  in  the  spatial  vicinity  of  the  earth's  surface,  so  that  the 
quantities  (2-34)  may  indeed  be  used  as  spatial  curvilinear  coor- 
dinates. 

The  matrix  M introduced  by  (2-19)  and  (2-20)  may  be 
written  as 


(2-35) 


it  is,  therefore,  nothing  else  than  the  Jacobian  matrix  of  the 
transformation  (2-32).  It  is  well  known  that  the  inverse  matrix 
M ^ is  then  simply  the  Jacobian  matrix  of  the  inverse  trans- 
formation (2-33)  : 


(2-36) 


This  may  also  be  shown  directly:  we  have 


3Yi  3Xj  _ 3y^ 


(2-37) 
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by  the  chain  rule  of  differential  calculus;  furthermore 

; 1 if  i = k , 

0 if  i ^ k ; 


- 6 
- ik 


(2-38) 


for  instance,  clearly 


3y. 
9y  ^ 


= 1 


9y. 

W- 


Therefore,  (2-37)  becomes 


3Yi  ax. 
ax.  ay^ 


6 . 


ik 


(2-39) 


which,  by  (2-35)  and  (2-36),  is  nothing  but  the  equation 


MW  = I 


(2-40) 


in  index  notation,  _I^  denoting  the  unit  matrix. 

Now  the  vector  m , defined  by  (2-29),  becomes  in  index 
notati on : 


ax. 

^ "a  ’ 


(2-41) 


9Xi  3Yj 

9T 

9^  "j  ’ 


(2-42) 


again  be  the  chain  rule.  Hence  (2-30)  becomes 


T 9 T j. 

" - ''1  n-  = f 


where  we  have  used  the  abbreviation 


(2-43) 


f = aW  + m 


(2-44) 


An  even  greater  simplification  is  achieved  by  introducing 
"quasi -spherical  coordinates"  p,  ({) , X by 


- ^ COS(|)  cosx 
p 

- coS(()  si  nx 

P 

- \ sin*  . 


(2-45) 


because  the  vector 


coord i nate 


Here  * and  x are  normal  latitude  and  longitude  as  before, 
because  the  vector  is  nothing  else  than  normal  gravity.  The 

coordinate  p is  taken  as  positive.  If  the  reference  ellipsoid 
becomes  a sphere,  then  p becomes  proportional  to  the  radius 
vector,  as  we  shall  see  below,  so  that  p,  <j) , x become  spherical 
coordinates;  hence  the  name,  quasi -spheri cal  coordinates. 
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Now 

il  = U- 

3p  8y^  3p 


again  by  the  chain  rule; 


9Yi  2 2 

T — = —T  cos<t>  cosx  = Y, 

9 p p 3 p i 

and,  generally, 

!li  = - 2 

9p  P 

by  (2-45),  Thus  (2-46)  becomes 
ll  = - 2 aJL 

9p  p ^i9Yj^  ’ 


(2-46) 


(2-47) 


and  (2-43)  reduces  to 

o|l  + 2T  = 2f  (2-48) 

It  should  be  pointed  that  (2-48),  in  spite  of  its  simplicity, 
is  rigorously  equivalent  to  (2-30);  there  is  no  further  approximation 
i nvol ved . 

What  is  the  geometrical  meaning  of  the  derivative  9T/9p  ? 

According  to  the  definition  of  a partial  derivative,  3/9p  means 
differentiation  with  respect  to  one  coordinate  p , the  two  other 
coordinates  <(> , \ being  held  constant.  This  means  differentiation 


i 

i 

1 

1 
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lines,  then  the  isozenithals  would  coincide  with  the  plumb  lines; 
as  the  normal  plumb  line  curvature  is  quite  small,  isozenithals 
and  plumb  lines  are  not  very  different. 

In  view  of  the  fundamental  importance  of  our  boundary 
condition,  let  us  approach  it  from  still  another  angle.  Let  t de- 
note the  arc  length  of  the  isozenithal  line,  measured,  e.g.,  from 
the  ellipsoid  positive  upwards  (so  that  it  represents  the  height 
above  the  ellipsoid  , measured  along  the  isozenithal).  Then  3/3t 
represents  a derivative  along  the  isozenithal,  in  the  same  way  as 
3/3p  . Therefore,  these  two  derivatives,  having  the  same  direction,, 

can  only  differ  in  scale,  that  is,  they  must  be  proportional: 


= 

3p 


(2-53) 


c = lx  : ix 

3t  9p 


1 1 lY 

"7^7^ 


(2-54) 


Hence,  by  (2-50) 


- 2(-  lX)"ii_ 


(2-55) 


and  the  boundary  condition  (2-48)  takes  the  form 


aT  1 9y  T - 

9 T Y 3 T 


1 lx  f 

Y 9t  ^ 


(2-56) 


The  right-hand  side  may  be  transformed  as  follows.  By  (2-44) 


we  have 


f = aW  + 


(2-57) 


Let  us  have  a closer  look  at  the  vector  m 
To  this  effect,  let 


X = x(t) 


(2-58) 


be  the  equation  of  the  isozenithal.  Then  the  vector 


X - -JT 


(2-59) 


will  be  the  unit  tangent  vector  of  this  curve  (it  will  be  a unit 
vector  since  t is  the  arc  length).  Then 


(2-60) 


e’’gnad  1 i± 

— ^ 3X . Or  3t 

1 


by  the  chain  rule.  Hence  there  follows  from  (2-42),  (2-47),  (2-55) 
and  (2-60); 


m^^grad  T = - y . 


3y.  'i 


1 3T 
7^3^ 


(1  9y\-13T 

S 9^^  3t 


= - (7  T 


(2-61) 


Since  the  vector  grad  T can  have  any  direction,  there  must  be 


= - (7  • (2-62) 

Hence  the  vector  m is  tangent  to  the  isozenithal;  since  t is 
positive  upwards,  the  negative  sign  implies  that  m is  directed 
downwa  rds . 

Thus 


= -(7 


is  nothing  else  than  the  component  of  the  gravity  vector  A£  in  the 
downward  direction  of  the  isozenithal.  Since  this  direction  is  very 
nearly  vertical,  Ag'  is  almost  equal  to  the  usual  gravity  anomaly 
Ag  in  the  sense  of  Molodensky. 

In  view  of  (2-63)  and  (2-64),  eq.  (2-57)  becomes 


f ■ 4“  M:) 


and  (2-56)  may  be  written  as 


il  - i ft  T = - ig' 
9 T Y 3 ^ 


(2-65) 


_1  9 Y 
Y 9 T 


(2-66) 


This  form  of  the  basic  boundary  condition  is  rigorously 
equivalent  to  the  preceding  forms  (2-30),  (2-43)  and  (2-48).  Though 
it  looks  less  simple,  it  is  very  important  because  it  allows  a 
comparison  with  the  form  in  which  the  boundary  condition  for 
Molodensky's  problem  was  usually  presented  earlier.  Take,  for 
instance,  eq.  (8-24b)  of  (Heiskanen  and  Moritz,  1967  , p. 300): 


9T  1 3y  t 

FIT  " Y 


- Ag 


(2-67) 


Here  the  derivative  9/9h  is  taken  along  the  normal  plumb  line. 

This  equation  involves  certain  approximations  (cf.  ibid. , p . 8 5 ) , 
which  are  practically  permissible  but  theoretically  not  rigorous. 

It  was  the  merit  of  T.Krarup  to  have  shown  that  (2-67)  becomes 
theoretically  exact  if  the  direction  of  the  normal  plumb  line  is 
replaced  by  the  direction  of  the  normal  isozenithal  (the  second  term 
on  the  right-hand  side  of  (2-66)  vanishes  if  the  telluroid  is 


defined  by  U 


as  usual ) . 
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The  boundary  condition  (2-66)  is  valid  on  the  telluroid  I , 
which  is  a known  surface.  The  problem  is  to  solve  Laplace's  equation, 
aT  = 0 , outside  I with  the  boundary  condition  (2-66).  Since  the 

isozenithal  is,  in  general,  not  normal  to  the  surface  I , we  have 
an  oblique  derivative  problem.  Such  problems  are  considerably  more 
difficult  than  boundary-val ue  problems  involving  normal  derivatives, 
such  as  Stokes'  problem. 

Spherical  Approximation.  - If  the  reference  ellipsoid  is  a 
nonrotating  sphere,  then 


Y 


(2-68) 


where  G is  the  gravitational  constant,  M the  total  mass,  and  r 
the  radius  vector  from  the  center  of  the  sphere  to  the  point  under 
consideration.  The  normal  gravity  vector  is  then  given  by 


X = ■ 


where 


e 


COS(j>  cosx 
coscfi  sinx 
s i n,), 


(2-69) 


(2-70) 


denotes  the  unit  vector  in  the  direction  of  the  radius  vector,  (|) 
and  X being  geocentric  latitude  and  longitude.  The  quantities  r, 
^ , X are  the  usual  spherical  coordinates. 

The  cartesian  components  of  j '"^y  thus  be  written 
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GM 

— COScji  COSX 


GM 

— cos<))  sinx  , 
r" 


(2-71) 


GM  . ^ 
sin(|> 

r 


The  comparison  with  (2-45)  shows  that  now 


= r//^  . 


(2-72) 


so  that  p is  r apart  from  a scale  factor. 

For  the  non-rotating  sphere,  the  plumb  lines,  as  well  as  the 
i sozen i tha 1 s , coincide  with  the  spherical  radii.  Thus,  now 


3 _ 9 

3t  " 3r 


(2-73) 


and 


A 9 Y = A 9 _ A 

Y 3t  y 3*"  ” 


by  (2-53).  Hence  (2-66)  reduces  to 


(2-74) 


|I  . 2 T . 
3 r r 


- Ag  + - aW 


(2-75) 


equivalent  to  (2-48)  but  with  the  right-hand  side  given  explicitly. 

The  boundary-va 1 ue  problem  expressed  by  Laplace's  equation 


at  = 0 


(2-76) 


29 


and  the  boundary  condition  (2-75)  in  spherical  coordinates  has  been 
called  by  Krarup  the  simple  Molodensky  problem;  it  is  the  one 
considered  in  virtually  all  practical  solutions  of  the  geodetic 
boundary  value  problem. 

The  reason  is  that,  although  the  reference  ellipsoid  is  not 
exactly  a sphere,  its  flattening  is  very  small,  about  0.3  %,  so  that 
on  tolerating  an  error  of  this  order  of  magnitude  in  equations 
relating  quantities  of  the  anomalous  gravity  field,  for  instance, 
in  the  boundary  condition,  we  can  formally  use  spherical  boundary 
condition  even  in  the  geodetic  case  of  a reference  ellipsoid.  This 
is  the  so-called  spherical  approximation;  for  a more  detailed 
explanation  cf.  (Heiskanen  and  Moritz,  1967 , pp . 87-88) . 

As  Krarup  has  pointed  out,  the  spherical  approximation  may 
be  interpreted  geometrically  as  the  mapping  of  the  actual  point  P 
into  an  auxiliary  point  P'  by  relating  the  quas i -spher i ca 1 co- 
ordinates p,  (() , X of  P to  the  spherical  coordinates  r',  ((i ' , x' 
of  P'  by 

r’  = , ,j)‘  = <1)  , X'  = X . (2-77) 

This  mapping  would  even  be  rigorous  if  also  the  Laplace  equation  were 
transformed  appropriately;  the  approximation  amounts  to  the  use  of 
the  untransformed  Laplace  equation.  The  reader  may  find  it  interesting 
to  compare  this  with  the  ideas  to  be  presented  in  sec.  8. 

As  we  have  already  said,  the  simple  boundary  condition  (2-75) 
has  been  used  in  almost  all  practical  solutions  of  the  geodetic 
boundary  value  problem.  This  is  already  true  for  Stokes  problem,  the 
gravimetric  determination  of  the  geo  id.  In  fact,  for  the  reference 
ellipsoid  itself,  which  is,  by  the  spherical  approximation,  mapped 
onto  the  sphere  r = R , we  have  with  aW  = 0 : 


L 


(2-78) 


I 
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which  is  the  boundary  condition  for  Stokes'  problem  (Heiskanen  and 
Mori tz , 1967 , p . 88) . The  solution  is  given  by  the  well-known  Stokes 
i ntegral . 

But  also  almost  all  practical  solutions  of  Molodensky's 
problem  presented  and  applied  so  far  are  based  on  the  spherical 
approximation,  beginning  with  (Molodenskii  et  a 1 . , 196  2 , pp . 1 18- 124 ) ; 
solutions  by  Arnold,  Brovar,  Marych,  Moritz,  Pellinen  and  others. 
For  a review  of  them  see  (Moritz,  1966  and  1969).  The  ellipticity 
has  been  taken  into  account  in  work  by  Zagrebin,  Molodensky, 
Bjerhammar  and  Koch;  Lelgemann  has  shown  that  the  effect  of 
ellipticity  on  geoidal  heights  and  deflections  of  the  vertical  is, 
in  fact , very  smal 1 . 

Let  us,  finally,  mention  that  there  are  two  senses  in  which 
"linear"  and  "nonlinear"  are  used  in  connection  with  Molodensky's 
problem.  In  theoretical  work,  such  as  the  present  report,  "linear- 
ization" is  with  respect  to  the  disturbing  potential  T or  the 
height  anomaly  ^ ; the  neglected  quantities  are  on  the  order  of 


10"^° 


which  is  certainly  always  negligible,  the  present  accuracy  being 
higher  than  10"^  only  in  rare  instances. 

Hence,  for  practical  applications,  the  "linear"  Molodensky 
problem  in  the  present  sense  is  practically  always  sufficient.  Higher 
approximations  and  their  convergence  are,  however,  of  basic  import- 
ance for  a rigorous  mathematical  investigation  of  the  existence 
and  the  uniqueness  of  a solution  to  Molodensky's  problem. 

The  other  sense  of  "linear"  and  "nonlinear"  is  used  in  more 
practically  oriented  work,  with  respect  to  terrain  inclination 
tan  e , because  the  practical  solutions  to  Molodensky's  problem 
(it  is  usually  the  "simple"  problem  mentioned  above)  are  based  on 


a series  expansion  with  respect  to  tan  6 . (It  is  in  this  sense 

that  "linear"  and  "nonlinear"  are  used  in  (Moritz, 1966  and  1969).) 
Since  tan  6 may  reach  rather  large  values,  "linear  solutions" 
in  the  second  sense  may  be  no  longer  adequate  in  mountainous  areas. 


3 . Hbrmander's  Linearization 

As  a preparation  to  his  fundamental  study  of  the  nonlinear 
Molodensky  problem,  Hbrmander  (1975)  first  transcribes  Krarup's 
linearization  process  into  modern  symbolic  notation. 

The  purpose  of  this  notation  is  to  exploit,  as  much  as 
possible,  the  analogy  between  linear  operators  and  ordinary  linear 
functions,  between  nonlinear  operators  and  ordinary  nonlinear 
functions,  etc.  In  this  way  it  is  possible,  for  instance,  to  develop 
differential  calculus  in  a unified  way  equally  valid  for  functions 
of  one  variable,  for  functions  in  n-dimensional  space  and  for 
functions  in  a Hilbert  or  Banach  space  (which  are  nothing  else 
than  linear  and  nonlinear  operators).  One  can  also  give  implicit 
and  inverse  function  theorems  that  are  equally  valid  for  functions 
of  one  or  several  variables  and  for  nonlinear  operators.  This  explains 
the  importance  of  such  a modern  notation  for  the  geodetic  boundary- 
value  problem. 

The  reader  will  find  presentations  of  this  "modern  analysis" 
in  books  such  as  ( D i eudonn§ , 1960)  and  (Loomis  and  Sternberg , 1968 ) , 
the  first  being  very  clear  but  rather  abstract,  the  second  more 
intuitive  and  accessible. 

In  this  notation,  vectors  and  functions  are  not  distinguished 
by  notation  from  ordinary  numbers.  Vectors  will  be  written  simply 
a or  b instead  of  £,  b (or  of  a^,  b^  in  indices  notation); 
functions  are  denoted  simply  by  f or  g instead  of  f(x)  or 
g(u,v)  . Thus  structural  similarities  between  numbers,  vectors 
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and  functions  are  stressed,  but  it  must  be  specified  at  each  instance 
what  these  symbols  really  denote. 

The  following  standard  mathematical  notations  will  also  be 

used : 

R ..  line  of  real  numbers  x , - «><x<“  ; 

R^...  three-dimensional  Euclidean  space. 

We  also  denote  by  a the  unit  sphere  in  R^  (we  do  not 

2 

use  the  standard  notation,  S , to  avoid  confusion  with  the  earth's 
surface  S ) . There  i s 

CT  = { xfR^  ; + X2  + Xj  = 1}  , (3-1) 

which  means  that  the  unit  sphere  is  the  set  of  all  points  of  R^ 
for  which  x^  + x^  + x^  = 1 holds. 

The  potential  W is  a function 

W : R^  ^ R , (3-2) 

which  means  that  the  function  W(x^,  X2 , X3 ) associates  to  each 
point  (x^,  X2 , X3 ) of  R^  a number,  namely  the  value  of  the 
gravity  potential  W at  this  point,  and  this  number  is  an  element 
of  R . In  this  sense,  the  function  W does  indeed  map  R into 
R , mapping  into  R being  understood  as  associating  a numerical 
value  (a  real  number). 

The  vector  (x^,  X2 , X3)  will  be  denoted  by  x , so  that 
we  may  also  write 

W = W(x) 

= V(x)  + ^^(Xj  + X2)  , 


(3-3) 
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where  V is  the  gravitational  potential  and  the  last  term  the 
potential  of  centrifugal  force,  u being  the  angular  velocity  of 
the  earth's  rotation. 

The  gravity  vector  g (no  underlining!)  is  the  gradient 

of  W : 


g = W 


(3-4) 


where  the  prime  denotes  differentiation;  since  x means  (x^, 
x^)  , it  is  quite  natural  that  W designs  the  gradient: 


8W  3W  3W 
3X^  ’ 3X2  ’ 3X2 


grad  W 


(3-5) 


How  can  the  physical  earth's  surface  S be  defined  in  the 
modern  way?  It  is  a mapping 


S : a ^ 


(3-6) 


What  does  this  mean?  Nothing  else  than  the  usual  parametric  re- 
resentation  of  S in  terms  of  the  parameters  $ (astronomical 
latitude)  and  A (astronomical  longitude): 

~ 3^($,A)  , 

X2  = S2($,A)  . (3-7) 

x^  = S2($,a)  , 

where  S^,  S^,  are  certain  functions  of  4>  and  A , or  briefly 


X = S ('t  ,A  ) 


(3-8) 
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I 
I 

j 

j 

In  fact,  * and  A define  a point  on  the  unit  sphere  a (as  we 

know  well  from  spherical  astronomy!),  and  x^,  x^,  x^  determine  , 

a point  in  . We  are  purposely  using  the  same  symbol  S for  the  i 

surface  and  its  equation  (3-7)  or  (3-8);  the  usefulness  of  this 

notation  will  be  seen  later. 

If  the  functions  (3-7)  are  differentiable  ( as  will  be 
assumed),  then  the  mapping  (3-6)  will  be  differentiable:  the  surface 
S is  a "differentiable  embedding"  of  the  unit  sphere  in  . 

On  putting 

f = u^  , A = U2  (3-9) 

we  can  write  (3-8)  in  the  form 

X = S(u^,  u^)  (3-10) 

or  still  more  briefly 

X = S(u)  , (3-11) 

where 

u = (u^,  u^)  = ($,A)  . (3-12) 

In  sec.  1 we  have  denoted  the  gravity  potential  of  S (in 
modern  terminology,  the  restriction  of  W to  S ) by  H . 

Clearly,  ¥ is  a mapping 

¥ : a ^ R , (3-13) 


which  simply  means  that  W , as  a function  defined  on  a surface,  is 
a ( real -val ued ) function  of  the  two  surface  parameters  $ and  A 


defining  a point  on  the  unit  sphere  a 

In  the  same  way,  the  restriction  'g  of  the  gravity  vector 
g to  S is  a function,  or  mapping, 

g ; a ^ (3-14) 

(R^  because  of  the  three  components  of  "g  ). 

How  can  we  relate  the  surface  restrictions  ¥ and  "g  to 
their  spatial  counterparts  W and  g ? In  usual  notation  we  have 

¥($,A)  = W(S^($,A),  S2($,A),  S3($,a))  . (3-15) 

This  means  that  this  restriction  to  F is  obtained  by  substituting 
the  surface  equation  (3-7)  into  the  spatial  function  W(xi,  X2 , X3) 
More  briefly  this  is  written 

W = W(S(u))  . (3-16) 

It  is  thus  a composite  function  of  u , that  is,  of  4>  and  A 

The  modern  notation  for  a composite  function  or  "function 
of  a function" 

f(^(x))  (3-17a) 

1 s 

fo(j)(x)  or  fo(|)  (3-17b) 

(read  "o"  as  "circle"  or  "composed  by"),  so  that  (3-16)  may  be 
written 


3-18 


36 


or  briefly 

y = WoS  . (3-19) 

In  fact,  S is  a mapping  a i by  (3-6),  and  W is  a mapping 

-V  R , by  (3-2),  so  that  the  composition  of  the  two  mappings 
(read  the  right-hand  side  of  (3-19)  from  the  right  to  the  left!)  is 
a mapping  o ^ R , as  it  should  be  by  (3-13). 

Now  it  is  clear  that  similarly 

g = goS  = W'oS  (3-20) 

(bj  3-4).  In  other  terms,  restriction  to  the  surface  S is 
equivalent  to  composition  by  the  functional  symbol  S . This  shows 
the  usefulness  of  the  present  notation,  using  S in  both  meanings. 

Linearization.-  To  proceed  with  the  linearization,  we  consider 
S,  W,  g,  y and  g as  smooth  (that  is,  sufficiently  often  differ- 
entiable) functions  of  a parameter  6 . Let  6=0  correspond  to 

the  normal  potential  and  to  the  telluroid,  and  e = 1 to  the 
actual  potential  and  to  the  physical  earth's  surface. 

In  other  terms,  let 

x=S(u;  0)  (3-21) 

denote  a set  of  surfaces  which  depend  smoothly  on  the  parameter  e , 

0 < e < 1 . The  limiting  surfaces  of  this  set  are: 

X = S(u;0)  ...  telluroid  \ , (3-22) 

X = S(u;l)  ...  earth's  surface  S ; 


for  0 < 0 < 1 we  get  intermediate  surfaces. 
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Similarly,  let 
W = W(x;0) 

denote  a set  of  gravity  potential  functions.  The  limiting 
of  this  set  are: 

W = W(x;0)  ...  normal  potential  U ; 

W = W(x;l)  ...  actual  potential  W 

Finally,  let 

¥ = ¥(u;e) 

denote  the  restriction  of  W(x;9)  to  the  surface  S(u;e) 

¥(u;e)  = W(S(u;0)  ;0) 

The  limiting  functions  in  this  case  are 
¥ = ¥(u;0)  , 

which  denotes  the  normal  gravity  potential  U on  the  telluroid,  and 
¥ = ¥(u;l)  , 

which  denotes  the  actual  gravity  potential  W on  the  physical 
earth's  surface. 

Similarly  we  proceed  with  the  gravity  vector 


(3-23) 

functions 

(3-24) 


(3-25) 


(3-26) 


g = g(u;0)  ; 


(3-27) 
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in  this  case 

g(u;0) 

denotes  the  normal  gravity  vector  x telluroid  and 

?(u;  1) 

I 

denotes  the  actual  gravity  vector  £ on  the  earth's  surface. 

The  differentiation  with  respect  to  0 will  be  denoted  by 

a dot: 

: F = If  . (3-28) 

i 

! To  establish  the  relation  to  the  preceding  sections  we  note 

that  AS  in  sec.l  corresponds  to 

1 

I AS  = S(u;l)  - S(u;0) 

i = Sab  , 

according  to  Taylor's  theorem.  However, 

A0  = 1 - 0 = 1 , (3-29) 

so  that  we  get 

AS  = S = £ . (3-30) 

In  other  words,  S is  nothing  else  than  AS  as  used  in  sec.l;  it 
also  equals  the  vector  £ given  by  (2-15)  since  for  fixed  u , 

S(u;l)  represents,  by  (3-21),  the  position  vector  x of  point  P , 

L _ J 
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and  S(u;0)  that  of  Q (Fig. 2). 

Similarly  we  get  for  the  (vector)  gravity  anomaly 

^■2-  “ 5.  earth's  surface  " ^ telluroid 
= g(u;l)  - g(u;0) 

= gA0 


or 


A£  = g , (3-31) 

so  that  our  former  A£  is  nothing  else  than  "g  in  the  present 
notation. 

Thus,  linearization  is  now  equivalent  to  differentiation 
with  respect  to  the  parameter  6 

Let  us  first  differentiate  ¥ as  given  by  (3-26): 

¥(u;0)  = W(S(u;0);0)  , (3-32) 

obtained  by  substituting  the  surface  equation 

X = S(u;e ) (3-33) 

into  the  spatial  expression 

W = W(x;0)  . (3-34) 

Since  0 enters  into  the  right-hand  side  of  (3-32)  in  two  ways, 
directly  and  also  indirectly  through  S(u;0)  , we  must  apply  the 

chain  rule: 
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¥(u  ;0  ) = W(S(u;0 ) ;0  ) + 

+ <grad  W(S(u;0);0)  , S(u;0)>  . (3-35) 

What  does  this  mean?  The  first  term  on  the  right-hand  side  denotes 
the  derivative,  with  respect  to  0 , of  the  function  W(x;0)  , 

disregarding  the  dependence  of  x = S(u;0)  on  0 . The  second 

term  takes  into  account  precisely  this  latter  dependence.  We  thus 
have  to  differentiate  W with  respect  to  x , getting 

W'  = grad  W (3-36) 

and  then  x = S(u;0)  with  respect  to  0 , obtaining 

X = S(u;0)  . (3-37) 

Now  both  grad  W and  (3-37)  are  vectors  in  , and  < , > is  to 

denote  the  inner  product  of  two  such  vectors: 


<a,b>  - a^b^  + a^b^  + a^b^ 

= a^b^  (3-38) 

in  index  notation.  In  fact,  the  last  term  in  (3-35)  might  also  be 
written  in  the  more  familiar  index  notation 


9W  • ^ 9^  ^ 

9X^  ” 9X^  90 


(3-39) 


which  expresses  the  usual  chain  rule. 

After  having  clarified  the  meaning  of  (3-35),  we  omit  all 
arguments  but  are  careful  to  denote  restriction  to  the  surface  S 
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¥ = iios  + <w'oS,  s> 

• • 

= WoS  + <goS,  S> 

• • 

= WoS  + <g,  S> 

In  an  analogous  manner  we  get 
g = goS  + <g  oS,  S> 

What  does  g'  mean?  We  have 


g 

g’ 


W'  = grad  W = 
W"  = 


3W 

3x. 

1 


3^W 

3 X . 3 X . 
1 j 


9 


(3-40) 

(3-41) 

(3-42) 

(3-43) 


which  is  nothing  else  than  the  second-order  gradient  tensor.  For 
6 = 0,  W reduces  to  the  normal  potential  U , so  that  then  g'  is 
nothing  else  than  the  matrix  M given  by  (2-19),  and  g'oS  is 
this  matrix  M taken  at  the  surface  of  the  telluroid.  Note  that 
S here  denotes  S(u;e)  ; that  is,  any  surface  between  the  telluroid 
(e  = 0)  and  the  earth's  surface  (6=1)  , including  these  two 

surfaces . 

We  must  now  find  out  the  meaning  of  W and  ¥ . Using 

(3-29)  we  have 

W = WA9  = W(x;  1)  - W(x;0) 

= W - U = T , 

so  that 


W = T 


(3-44) 
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and  WoS  means  the  disturbing  potential  T at  the  surface  of  the 
tel  1 uroid . 

« 

On  the  other  hand,  ¥ means  something  different.  We  have 

¥=yA0=¥(u;l)-¥(u;O)  , 

which  is  the  difference  between  the  actual  potential  W at  the 
, earth's  surface  and  the  normal  potential  U at  the  telluroid,  in 

other  words,  using  the  notation  of  sec. 2,  the  difference  W^  - 
I which  by  (2-3)  is  aW  . Thus 

f 

j 

I ¥ = aW  . (3-45) 

* 

I We  finally  have 


g = gAe  = g(x;l)  -g(x,0) 

= grad  W - grad  U = grad  T . (3-46) 

Let  us  now  collect  equations  (3-30),  (3-31),  (3-44),  (3-45) 
and  (3-46): 

S = 1 , 

W = T , 

g = grad  T , (3-47) 

¥ = AW  , 
g = A^  . 


These  equations  form,  so  to  speak,  part  of  a dictionary 
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notation  and  vice  versa. 

We  thus  wee  that  (3-40)  is  the  "translation"  of  (2-25), 
whereas  (3-41)  corresponds  to  (2-26). 

We  now  proceed  as  in  sec. 2.  We  solve  (3-41)  for  S : 

S = (g'°s)-^t  - goS)  . (3-48) 

assuming  that  the  "Marussi  condition", 

det(g 'oS)  i 0 

is  satisfied,  that  is,  that  the  3x3  matrix 

M = g'oS  (3-49) 

is  invertible.  This  is  always  the  case  for  e = 0 , for  the  normal 

ellipsoidal  potential  U 

The  substitution  of  (3-48)  into  (3-40)  then  gives  the 
desired  boundary  condition.  We  first  obtain 

W = WoS  + <goS, (g'oS)"^ (g  - goS)>  . (3-50) 

Generally  we  have  for  3-vectors  a and  b and  for  a 3x3  symmetric 
matrix  C : 

<a,Cb>  = a'^Cb  = b'^Ca  = 

» 

= <b,Ca>  = <Ca,b>  , (3-51)  i 

) 

using  matrix  notation  in  between  and  taking  into  account  the 
symmetry  of  the  matrix  C and  also  the  symmetry  of  the  inner 
product  < , > 
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By  means  of  this  formula,  eq.  (3-50)  becomes 

^=WoS+<(g'oS)^goS,g-goS>  , (3-52) 

We  now  put 

m = - (g')"^g  (3-53) 

This  3-vector  m has  already  been  used  in  sec. 2;  cf.  (2-29);  it  is 
directed  along  the  isozenithal.  Thus  (3-52)  becomes 

* * * 

= WoS  - <moS,g-goS> 

• ^ . 

= WoS  - <moS,g>  + <m,g>oS  ; (3-54) 

as  oS  means  restriction  to  S , there  is 

<moS,goS>  = <m,g>oS  . (3-55) 

On  rearranging  we  get 

(W  + <:m,g>)oS  = ¥ + <moS,'g>  . (3-56) 

Let  us  differentiate 

g(x;9)  = grad  W(x;9) 

with  respect  to  e . Since  the  differentiations  with  respect  to 
0 and  to  X (gradient)  are  independent  and  can  be  interchanged,  we 
get 


g = grad  W 


(3-57) 
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Thus  (3-56)  becomes 

(W  + <m,gracl  W>)oS  = "W  + <nioS,'g>  . (3-58) 

This  condition  differs  from  (2-30)  only  by  the  notation;  cf.,  in 
particular,  the  "dictionary"  (3-47). 

It  is  appropriate  to  take  the  centrifugal  potential  in 
space  as  independent  of  the  parameter  e . Then  only  the  purely 
gravitational  potential  V depends  on  e , and  W depends  on  e 
only  through  V . Thus 

W = V (3-59) 

and  (3-58)  takes  the  final  form 

(V  + <m,gradV>)oS  = if  + <moS,'g>  . (3-60) 

This  is  the  fundamental  boundary  condition  in  the  form  given  to  it 
by  Hbrmander. 

• • 

The  form  (3-60 ) clearly  brings  out  the  fact  that  T = V = U 
is  a ^'^.^monic  function,  satisfying  Laplace's  equation 

aV  = 0 . (3-61) 

This  is  less  obvious  in  (3-58)  since  V but  not  W is  a harmonic 
function . 

It  is  very  interesting  to  compare  the  equivalent  forms 
(2-30)  and  (3-60).  There  is  no  doubt  that  (2-30),  its  notations 
and  the  operations  leading  to  it  are  much  more  familiar  to  the 
geodetic  reader:  we  all  know  what  T and  Ag  are  and  are  all 
familiar  with  vector  and  matrix  operations. 


w 
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Equation  (3-60)  and  the  operations  leading  to  it  are  much 
more  abstract,  apart  from  being  unfamiliar  to  most  geodesists  at  the 
present  time  (will  this  change  in  the  near  future?).  On  the  other 
hand,  the  new  notation  possesses  considerable  logical  consistency  and 
rigor.  We  need  less  special  symbols  (no  T , no  c etc.),  restriction 
to  a surface  is  clearly  expressed,  and  differentiation  with  respect 
to  9 replaces  A and  other  differences.  It  is  sufficient  to  look 
at  (3-47);  note  how  clearly  the  new  notation  distinguishes  between 
T = W - U and  aW  = W - U ! 

P P P Q 

Furthermore,  (3-60)  is  formulated  in  terms  of  derivatives 
with  respect  to  e , rather  than  in  terms  of  approximate  finite 
differences.  In  the  new  notation  it  is  obvious  that  we  have  1 i near- 
ization  without  approximation;  this  is  less  evident  in  the  usual 
notation.  This  fact  makes  it  possible  to  linearize  not  only  at 
0=0  (which  would  be  the  first  stage),  but  also  at  any  intermediate 
0,  0 < 0 < 1 . (To  repeat,  S in  (3-60)  and  similar  equations 

denotes  S(u;9)  for  any  9 between  0 and  1 , including  9=0 

(telluroid)  and  0=1  (earth's  surface).)  If  we  linearize  at  9=0 
and  take  A0  = 1 , then  (3-60)  reduces,  in  fact,  to  (2-30),  but  it 

is  valid  for  linearization  at  any  0 and  is  thus  much  more  generally 
useful:  linearization  means  here  simply  differentiation. 

In  fact,  it  is  the  form  (3-60)  which  permits  the  method  of 
successive  approximation,  starting  at  0 = 0 and  taking  consecutively 
9 closer  and  closer  to  1 , which  will  be  outlined  in  sec. 5. 


4 . Existence  and  Uniqueness  of  Solution 
for  the  Linear  Molodensky  Problem 

We  shall  now  investigate  the  existence  and  uniqueness  of  the 
solution  of  the  linear  Molodensky  problem.  As  an  introduction  we 
examine  first  Stokes'  problem  and  the  "simple  Molodensky  problem". 
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following  (Mor i tz , 197 1 , pp . 27-3 1 ), before  considering  the  general 
linear  problem  following  ( Hbrmander , 1975 , cha pter  I). 

The  Problem  of  Stokes.-  Stokes'  problem  is  the  boundary- 
value  problem  in  its  simplest  form:  given  the  gravity  anomaly  on 
a sphere,  to  determine  the  anomalous  potential  T on  and  outside 
the  sphere,  assuming  T to  be  harmonic  outside  this  sphere.  The 
corresponding  ooundary  condition  is  (2-78);  since  the  radial  direction 
is  normal  to  the  bounding  sphere,  the  oblique-derivative  problem 
reduces  in  this  case  to  a problem  involving  normal  derivatives, 
which  is  much  simpler. 

The  solution  is  given  by  Stokes'  integral  formula 

T(e  ,A)  = Tq  + ^//AgS(^)da  + Ti  (0  ,X)  , (4-1) 

a 

which  expresses  T on  the  given  sphere  in  terms  of  Ag  on  this 
sphere.  Here  T^  is  a fixed  constant  related  to  the  mass  of  the 
earth,  and 

Tj  (0  ,X  ) = AjSinOcosX  + A2Sin0sinX  + f^^cosQ  (4-2) 

is  a spherical  surface  harmonic  of  the  first  degree.  Polar  distance 
0 and  longitude  X are  spherical  coordinates,  and  A^  , A^  , A^  are 
arbitrary  constants  which  have  the  following  physical  interpretation 
(Heiskanen  and  Mor i tz , 1967  , p . 99 ) . Let  €2*  denote  the 

rectangular  coordinates  of  the  earth's  center  of  gravity,  the  origin 
being  the  center  of  the  ellipsoid.  Then,  approximately, 

Ai=?5.  , (4-3) 

where  "g  denotes  a mean  value  of  gravity  over  the  earth.  Therefore, 
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nonzero  mean  that  the  center  of  the  reference  ellipsoid  does 

not  coincide  with  the  earth's  center  of  mass. 

A necessary  and  sufficient  condition  for  Stokes'  problem 
to  be  solvable  for  continuous  boundary  values  is  that  the  function 
Ag  does  not  contain  spherical  harmonics  of  the  first  degree.  In 
other  terms,  Ag  must  be  orthogonal  to  any  harmonic  function  of  the 
first  degree  Y^(0,X)  : 

//Ag(9 ,A)Y^(0 ,X)da  = 0 ; (4-4) 

a 

cf.  (Heiskanen  and  Mori tz , 1967  , p . 97  ) . Since  Yj(0,X)  contains 
three  constants,  this  equation  comprises,  in  fact,  three  independent 
conditions. 

The  solution  (4-1)  contains  three  free  constants  A^,  A^  . 

The  solution  can  be  made  unique  by  putting  all  A^  = 0 , which  means 

that  the  first-degree  harmonic  (4-2)  vanishes. 

The  fact  that  Ag  must  satisfy  three  conditions  and  that 
the  solution  (4-1)  contains  three  free  constants  expresses  the 
so-called  Fredholm  alternative;  see  below. 

It  should  also  be  pointed  out  that  a solution  (4-1)  with 
Ai  / 0 A2  , is  physically  impossible,  although  it  is  mathe- 
matically valid  as  a solution  of  the  boundary-value  problem  defined 
by  aT  = 0 outside  the  sphere  and  by  the  boundary  condition  (2-78) 
on  the  sphere. 

In  fact,  for  T to  be  harmonic  and  zero  at  infinity,  the 
centrifugal  potential  contained  in  both  W and  U must  be  equal, 
so  as  to  drop  out  in  T = W - U . This  requires  that  the  axis  of  the 
reference  ellipsoid  coincides  with  the  earth's  axis  of  rotation. 

If  this  common  axis  is  taken  as  x,  axis,  then  the  centrifugal 
potential  is 
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12,2^2.  / /I  c \ 

jw  (x^  + x^)  . (4-5) 

Indeed,  if  the  two  axes  were  only  parallel  and  separated  by  the 
vector 


(6Xj,6X2>0)  , 

then  T would  contain  a term 

a)^(x^6x^  + (4-6) 

due  to  the  difference  of  the  two  centrifugal  potentials;  this  term 
and  therefore  T , would  not  be  zero  at  infinity.  The  same  would 
hold  if  the  two  axes  were  not  parallel. 

So  the  two  rotation  axes  must  coincide.  Since  the  earth's 
rotation  axis  passes  through  the  center  of  mass  for  physical  reasons, 
and  since  the  axis  of  the  ellipsoid  contains  the  center  of  the 
ellipsoid  for  reasons  of  symmetry,  both  centers  must  lie  on  the 
common  axis,  which  is  taken  as  the  x^  coordinate  axis.  This  implies 
that  the  two  centers  can  differ  only  in  the  x^  coordinate,  so  that 
Cl  and  C2  > therefore  Ai  and  by  (4-3),  must  be  zero. 

Thus,  if  a solution  (4-1)  is  to  be  physically  meaningful, 
only  A^  can  differ  from  zero,  so  that  the  solution  for  a rotating 
earth  has,  in  reality^  only  one  degree  of  freedom.  Since  Aj  = A^  = 0 , 
it  is  quite  natural  to  take  also  A^  = 0 , thus  letting  the  center 

of  the  reference  ellipsoid  coincide  with  the  earth's  center  of  mass. 

The  Simple  Molodensky  Problem.-  This  is  the  linear  Molodensky 
problem  for  a spherical  reference  surface  (sec. 2).  We  shall  prove 
existence  and  uniqueness  of  the  solution  for  this  problem  by 
establishing  a one-to-one  correspondence  between  the  Stokes  problem 
and  the  simple  Molodensky  problem. 


variable  radius  vector  of  the  point  under  consideration  (ibid.  ,p.90) 
As  the  boundary  values  of  Ag  , and  hence  of  F , are  given  on 
the  surface  I , we  can  compute  F , and  hence  Ag  , at  every 
point  outside  ^ by  solving  an  external  Dirichlet  problem,  which 
is  uniquely  solvable  for  continuous  boundary  data  (cf.  Kel 1 ogg , 1929 , 
p.314).  In  particular,  this  gives  Ag  at  every  point  P between 
the  surfaces  I and  S'  --to  be  denoted  by  Ag(r)  --and  on  the 
Brillouin  sphere  itself--to  be  denoted  by  Ag' 

From  the  values  of  Ag  along  a radius  it  is  straightforward 
to  compute  radial  differences  of  the  potential  T : by  (Heiskanen 

and  Moritz,  1967,  p.92)  we  have 


l^(r^T)  = - r^Ag(r) 


which  on  integration  gives 


(r  ■r)p  - (r  T)p,  = /r  Ag(r)dr 

p 


or  with  the  symbols  of  Fig. 4, 


(4-8) 


(4-9) 


^pT  - RgT’  = /r^Ag(r)dr 


(4-10) 


T denoting  the  potential  on  the  telluroid  and  T'  on  the  Brillouin 
sphere . 

Now  we  can  solve  Molodensky's  problem  by  the  following  three 

steps : 

1.  Computation  of  Ag(r)  and  Ag'  by  solving  the  external 
Dirichlet  problem. 
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2.  Determination  of  T'  from  Ag'  by  solving  Stokes' 
problem  for  the  sphere  S' 

3.  Computation  of  T at  ^ from  T'  at  S'  by  (4-10). 

Steps  1 and  3 are  one-to-one  because  Dirichlet's  problem  is 
uniquely  solvable  and  because  different  functions  Ag  on  I corre- 
spond to  different  functions  Ag'  on  S'  and  vice  versa.  Thus  the 
question  of  solvability  of  Molodensky's  problem  for  the  telluroid 
1 has  been  reduced  to  the  question  of  solvability  of  Stokes'  problem 
for  the  sphere  S'  , to  which  the  answer  has  been  given  above.  For 
the  simple  Molodensky  problem,  therefore,  we  have  exactly  the  same 
situation  concerning  existence  and  uniqueness  of  solution  as  for 
Stokes'  problem:  Ag  must  satisfy  three  conditions,  which  may  be 
expressed  in  the  form,  analogous  to  (4-4), 

//Ag' (0,A)Y^(0,X)da  = 0 , (4-11) 

a 

which  means  that  the  upward  continuation  of  Ag  to  S'  must  not 
contain  any  first-degree  spherical  harmonic. 

Corresponding  to  these  three  conditions,  the  solution  for 
T'  , and  consequently  also  for  T , will  contain  three  free 
constants  (this  is  true  if  the  linear  boundary  value  problem  is 
considered  in  itself;  for  physical  reasons,  two  of  these  constants 
must  be  zero).  Again  we  get  a unique  solution  by  requiring  the 
spatial  function  T to  have  a form  that  contains  no  first-degree 
spherical  harmonics. 

The  Linear  Molodensky  Problem.-  The  general  linear 
Molodensky  problem  for  an  arbitrary  reference  surface  for  an  arbi- 
trary reference  potential,  as  formulated  in  sections  2 and  3,  is  an 
oblique-derivative  problem. 

The  classical  boundary  value  probl ems--the  Dirichlet  problem 
and  problems  involving  normal  deri vati ves--can  be  formulated  in  terms 
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of  Fredholm  integral  equations  of  the  second  kind,  and  the  well- 
known  Fredholm  alternative  holds  (cf.  Kel 1 ogg , 1929 , p . 298) : 

If  the  homogeneous  boundary-value  problem  has  no  non-zero 
solution,  then  the  corresponding  nonhomogeneous  problem  is  solvable 
for  arbitrary  continuous  boundary  values. 

If  the  homogeneous  problem  has  n independent  solutions, 
then  the  boundary  values  must  satisfy  n independent  conditions  for 
the  corresponding  nonhomogeneous  problem  to  be  solvable,  and  the 
solution  depends  on  n free  parameters  (because  of  the  n inde- 
pendent solutions  of  the  homogeneous  problem). 

An  example  is  furnished  by  Stokes  problem,  in  which  n = 3 . 

An  analogous  formulation  of  the  oblique-derivative  problem 
leads  to  singular  integral  equations  for  which  the  Fredholm  alter- 
native is,  in  general,  no  longer  valid.  An  example  is  Molodensky's 
integral  equation  which  is  no  longer  a Fredholm  equation  of  the 
second  kind  (contrary  to  what  is  sometimes  said  in  the  literature). 

However,  if  the  oblique-derivative  problem  is  regular,  that 
is,  if  the  direction  of  the  derivative  is  nowhere  tangential  to  the 
boundary  surface,  the  Fredholm  alternative  is  still  valid,  in  spite 
of  the  singularity  of  the  corresponding  integral  equation;  cf. 

(Mi randa , 1970 , p . 86) ; this  means  that  the  number  of  conditions  on  the 
boundary  data  f , given  by  (2-57),  is  equal  to  the  degree  of 
freedom  in  the  solution,  say  n 

In  the  simple  Molodensky  problem  we  again  had  n = 3 .In 
the  present  general  linear  case,  n must  be  at  least  three,  in  view 
of  the  three  degrees  of  freedom  in  the  spatial  shift  of  the  origin, 
but  perhaps  n = 4 or  5 ? 

Hdrmander  proved  that  even  in  the  general  form  of  the 
linar  Molodensky  oroblem  n equals  3 . First,  the  problem  defined 

by  (3-60)  and  (3-61)  is  reformulated  as  follows,  putting 
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V = T , (4-12) 

¥+<;moS,g’>  = f (4-13) 

(the  same  notations  were  used  in  sec. 2).  It  thus  hocomes:  to  determine 
a function  T satisfying 

1.  Harmonicity:  aT  = 0 outside  S . (4-14) 

2 . Boundary  condition  on  S ; 

(T  + <m,grad  T>)oS  = f . (4-15) 

3 . No  first-degree  harmonic; 

T(x)  = ^ + O(^)  , r ^ 00  , (4-16) 

r 

c being  some  constant. 

Hbrmander  proved  that  the  corresponding  homogeneous  problem,  that 
is,  (4-14),  (4-15)  and  (4-16)  with  f h 0 , has  the  unique  solution 

T = 0 . The  general  solution  of  the  homogeneous  problem,  without 

imposing  (4-16),  therefore  contains  the  three  independent  spherical 
harmonics  of  degree  1.  This  proves  that  n = 3 also  for  the  general 
linear  Molodensky  problem  (if  n were  >3  , the  solution  of 

(4-14),  (4-15),  (4-16)  would  no  longer  be  unique). 

The  principle  of  the  Hbrmander's  proof  of  uniqueness  may  be 
illustrated  by  the  corresponding  proof  for  Stokes'  problem.  Take 
Laplace's  operation  in  spherical  coordinates  r,  6,  X (cf.  Heiskanen 
and  Mori tz , 1967  , p . 19)  : 


at 


2 ar'  ir’ 
r 


1 


r s i ne  8 0 


^(sinefl)  . 


9^T 


~2  ; T~ 
r sin  0 3X 


(4-17) 
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F 

i 


Let  us  form  the  integral 

J = ///2r|^ATr^si nQdrdedA  (4-18) 

r > 1 

over  the  exterior  of  the  reference  surface,  which  we  identify  with 
the  unit  sphere  a by  choosing  the  radius  as  unit  of  length.  Since 
aT  = 0 outside  a , this  integral  will  be  zero: 

J = 0 . (4-19) 

Let  us  evaluate  the  integral  using  (4-17).  We  get 

’ (4-20) 

where  the  correspond  to  the  three  summands  on  the  right-hand 

side  of  (4-18).  We  shall  also  use  the  abbreviations 

dv  = r^sinedrdedx  (4-21) 

for  the  volume  element  and 

da  = sinededx  (4-22) 

for  the  surface  element  of  a , so  that 

dv  = r^drda  . (4-23) 


With  these  notations  we  have 


r 
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ffr(r^|I)  sinedrdedx 


r > 1 


There  is 


3r'’  dr' 


2rT  + r^T 
r rr 


Where  T denotes  dT/dr  as  usual,  and  further 


= 4r^T^  + r^|-(T2)  , 

r 8r'  r'  * 


SO  that 


+ /r^|-(T2)dr  . (4-25) 


1 1 ^• 


The  last  term  is  transformed  by  partial  integrati 


on : 


= (r^T^)_-  (.3^^,  JsrVdr 

1 * * J i 


(4-24) 


(4-26) 
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- 2 

The  first  term  on  the  right-hand  side  is  zero  since  = 0(r  ) , 

in  the  second  term  we  have  r = 1 , so  that  on  substituting  this 

expression  into  (4-25)  we  get 

CO  00 

/2rT  |-(r2T  )dr  = /r^T^dr  - (T^)  , . 

r3r  r'  r rr=l 

On  multiplying  this  equation  by  da  and  integrating  over  a , 
taking  (4-23)  into  account,  we  get  for  (4-24): 


J,  = /J/T^dv  - /jT^da  . (4-27) 

r>l  a 

Now  we  compute  . In  view  of  (4-17)  and  (4-18)  we  have 

'=2  ■ = 
r> 1 r s 1 ne 

= //j2rT^|^(sin0Tg)drdedX  (4-28) 

r > 1 

Let  us  first  perform  the  integration  with  respect  to  0 .By 
partial  integration  we  get 


/2rT^|^(sin0T0)d0 


= (2rT  sineT  ) 

' r 0 0 = Tr 

7T 

- (2rT  sin0T  ) - /2rT  sin0T  de 

' r 0 ' 0=0  ‘ r 0 0 


f 
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- JerT^gTgSinede 


/ryf:(TQ)sin0de  . 


Thus  (4-28)  becomes 


= - //Jr|7CrQ)drsin9dedx 


r>  1 


(4-29) 


The  integral  with  respect  to  r is  again  transformed  by  partial 
i ntegrati on : 


We  multi  ply  by  da  and  integrate  over  a , obtaining  in  view  of 
(4-23): 


r>  1 r a 


(4-30) 


The  integral  is  transformed  in  exactly  the  same  way 

as  , with  the  result 


J,  = lll^r^  Tjdv  * II 


9 J J J o 9 \ 

^ r>ir^sin^0  ^ 


. 2 
a s 1 n 0 


da 


(4-31) 
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and  we  are  through. 

We  only  have  to  collect  (4-E7),  (4-30)  and  (4-31)  and  to 
take  (4-19)  and  (4-20)  into  account.  The  result  is 


■2*2  Tj)dv  * 
r ^ s i n ^ e ^ 


+ //(Tq  + 


-4-  T")da  - //T"da 
sine  a 


(4-32) 


Now  is  nothing  else  than  the  radial  component  of  the  gradient 

grad  T , and 


It 

r e ~ ree  ’ 

(4-33) 

1 j ^ 3T 
rs 1 ne  A rs i ne  a A 


are  its  horizontal,  or  tangential,  components,  the  total  tangential 
component  having  the  squared  magnitude 


T 


2 

t 


2 • 

r s 1 n e 


(4-34) 


Thus  (4-32)  may  finally  be  written 


By  means  of  this  formula  we  can  prove  stability  for  Stokes' 
problem  as  follows.  The  boundary  condition  (2-78)  becomes  for  the 
homogeneous  problem  (Ag  = 0)  with  R = 1 : 

T + 2T  = 0 , (4-36) 

r 

from  which  we  get 

T = - 2T  on  CT  . (4-37) 

r 

If  T on  the  sphere  a is  expanded  into  a series  of  Laplace's 
surface  harmonics, 


T(0,A)  = lT^(0,x)  , 

O 


then 


= in(ntl)//T^d„ 
a o a 


(4-38) 


(4-39) 


cf.  (Molodenskii  et  al  . , 1962  , p . 87 ) . Furthermore,  by  (4-36), 

//T^da  = = 

a a 


= 4lffT^dG  , (4-40) 

° a ^ 


because  of  the  orthogonality  of  the  spherical  harmonics. 
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Since  the  first  term  in  (4-35)  can  never  be  negative,  there 

must  be 

//(T^  - T2)da  < 0 . (4-41) 

0 

This  equation  is  interesting  in  itself,  since  it  shows  that  the 
average  square  of  the  horizontal  gradient  will  always  be  smaller 
than  (or  equal  to,  but  only  for  T = 0 ) the  average  square  of  the 
vertical  gradient;  this  is  obviously  true  for  a sphere  of  any  radius 

R . 

Here  we  shall  use  (4-41)  for  the  proof  of  uniqueness.  We 
substitute  (4-39)  and  (4-40)  and  obtain 

00 

I [n(n  + l)  - 4]/|T^da  < 0 . (4-42) 

o a 


The  coefficients  within  brackets  are  positive  for  n = 2,3,4,  ... 
So  if  we  can  show  that  = 0 for  n = 0 and  n = 1 , then  all 
other  T^  must  also  be  zero  if  (4-42)  is  to  hold. 

For  n = 1 we  have,  in  fact,  = 0 because  of  (4-16). 
From  the  well-known  spheri cal -harmonic  relation  (cf , He i s kanen  and 
Mori tz , 1967  , p . 97 ) 


T 

n 


(4-43) 


we  get  for  n = 0 and  R = 1 


I 
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I 

\ 

i 

t- 


because  Ag  is  identically  zero,  for  the  homogeneous  Stokes  problem. 

Thus  we  have  proved  that 

T = 0 (4-44) 

is,  indeed,  the  only  solution  of  the  homogeneous  Stokes  problem  if 
(4-16)  is  prescribed. 

Of  course,  this  proof  is  considerably  more  involved  than  the 
simple  reasoning  concerning  Stokes  problem  using  spherical  harmonics, 
as  presented  at  tfie  beginning  of  this  section.  The  essential 
advantage  of  the  quadratic  condition  (4-35)  is,  however,  that  it 
can  be  generalized  to  the  linear  Molodensky  problem. 

This  has  been  done  by  Hdrmander.  His  proof  is,  however, 
extremely  involved  and  laborious  and  cannoc  be  given  here.  Even 
his  uniqueness  theorem  ( Hdrmander , 1975  , pp . 22-23)  is  so  complicated, 
involving  many  expressions  and  parameters,  that  it  cannot  be 
reproduced  in  the  present  report. 

Let  it  be  sufficient  to  mention  that  Hdrmander's  theorem 
contains  a number  of  parameters  which  depend  on  properties  of  the 
earth's  topography.  Larger  slopes  of  the  terrain  (say  60°)  are 
permitted  provided  they  do  not  occur  too  frequently.  Although  a 
detailed  study  of  fitting  Hdrmander's  parameters  to  the  actual 
earth's  topography  has  not  yet  been  made,  it  appears  that  the  theorem 
is  general  enough  to  ensure  the  uniqueness  of  solution  of  the  linear 
Molodensky  problem,  with  an  ellipsoidal  reference  field,  for  the 
actual  topography  of  the  earth,  the  exi stence  of  a solution  being 
generally  guaranteed  by  the  theory  of  the  oblique  derivative  oroblem. 


4 
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5 . Nash-Hormander  Iteration 

In  sec.l  we  have  outlined  Newton's  method,  which  is  widely 
used  in  advanced  implicit  and  inverse  function  theorems  of  nonlinear 
functional  analysis;  cf.  (Sternberg , 1969)  . 

In  his  treatment  of  the  nonlinear  Molodensky  problem, 
Hbrmander  (1975, chapter  II)  used  a different  iteration  method, 
which  is  a discrete  scheme  analogous  to  a continuous  method  used  by 
Nash  (1956). 

The  essence  of  the  two  methods  is  illustrated  by  Figures  5 
and  6.  The  problem  in  both  cases  is  the  same:  consider  a function 

y = $(x)  , (5-1) 

which  can  simply  be  a function  of  one  variable,  as  shown  in  the 
figures,  or  a nonlinear  operator  mapping,  for  instance,  one  Banach 
space  X into  another  Banach  space  Y , so  that 


$ : X > Y 


(5-2) 


i 

Given  the  value  't’(u)  of  $ at  some  unknown  point  x = u , to 

determine  u . ; 

We  proceed  as  follows,  writing  the  formulas  for  the  simple 
case  of  an  ordinary  function  of  one  variable;  the  Banach  space  case 
is  formally  quite  the  same.  We  assume  an  approximate  value  u^  for  ' 

u and  compute  ^(u^)  . The  further  procedure  is  different  for  the 

two  methods. 

In  Newton's  method  (Fig. 5)  we  intersect  the  tangent  at  | 

4>(u)withtheline  I 

o' 


y = $(u)  = const. 


(5-3) 
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parallel  to  the  x-axis;  this  line  is  known  because  4'(u)  is  given. 

In  this  way  we  obtain  u^  and  can  calculate  i})(u^)  . The  tangent 

at  ^(u^)  is  again  intersected  with  the  line  (5-3),  and  the  proce- 
dure is  repeated.  The  sequence  u , u^,  u^,  ...  tends  to  u provided 
the  iteration  converges. 

In  the  Nash-Hbrmander  scheme  (Fig. 6)  the  known  difference 


f = $ (u)  - $ (u^) 

is  split  up  in  a suitable  way  into  a convergent  series: 


(5-4) 


f = 

The  tangent 

y = 

which  gives 

y = 


f + f + f + 

o 1 2 


O 


at  $(u  ) is  now  intersected  with  the  line 

o 

$(u  ) + f = const., 
o'  o 

u^  . The  tangent  at  $(u^)  is  intersected 

<s>(u  ) + f + f,  = const., 

'o'  o 1 


(5-5) 


(5-6) 


with 


(5-7) 


which  gives  u^  . The  procedure  is  indefinitely  repeated. 

The  respective  iteration  procedures  are  illustrated  in  the 
figures  by  heavy  lines  with  arrows.  The  comparison  of  the  two 
figures  already  gives  an  indication  that  the  second  iteration 
scheme  may  have  certain  advantages  over  Newton's  method.  For  instance, 
in  the  Nash-Hbrmander  scheme,  the  u^  increase  more  or  less 
monotonical ly  to  u , whereas  in  Newton's  method  u^  may  be  "way 
out"  (Fig. 5).  In  the  second  scheme  we  have,  approximately. 


u 


k+  1 


- 1 


- u 


(5-8) 


so  that  the  convergence  of  the  approximation  can  be  controlled 
almost  arbitrarily  by  suitably  selecting  the  terms  in  (5-5);  in 
the  Newton  method,  the  convergence  is  quadratic  according  to  (1-22) 
The  formulas  for  Newton's  method  have  already  been  given  in 
sec.l.  In  the  sequel  we  shall  be  concerned  exclusively  with  the 
Nash-Hbrmander  iteration.  From  Fig. 7 we  read  off: 


Figure  7 

The  iteration  scheme 
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Here  we  have  put 

= [$'(\)]  '^  . (5-13) 

where  $'  denotes  the  derivative  with  respect  to  u . Geometrically, 
is  the  error  due  to  nonlinearity  of  the  function  $ ; it  is  the 

"height  of  the  steps"  in  Figures  6 and  7. 

According  to  Fig. 6 we  start  our  iteration  with 


h = f (5-14) 

o o 

by  (5-10),  since  there  is  no  error  ^ to  start  with.  Then,  by 
(5-11)  and  (5-12), 


6 

O 


= ip(u^)h 


O 


Next  we  put  k = 1 

e = ^(u,)  - $(u  ) 
o 1 ' o 


(5-15) 

(5-16) 


in  (5-9)  through  (5-12),  obtaining 
- h 

O 


h = f - e 
1 1 o 


6^  = i|<(u^)h^ 


U2  = u^  + S 


(5-17) 


The  procedure  in  repeated  with 


k = 2,3,.. 
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An  essential  feature  of  the  scheme  just  described  is  that 
the  error  committed  at  the  (k-l)th  step  is  taken  into 

account  and  corrected  at  the  next  step  by  means  of  (5-10). 

Unfortunately,  this  very  simple  iteration  scheme  cannot  be 
proved  to  converge  for  the  Molodensky  problem.  Therefore,  Nash  and 
Hbrmander  introduced  a suitable  smoothing  process,  which  also 
provides  a natural  determination  of  the  terms  f^^  in  (5-5). 

Let  us  i ntroduce  a sequence  of  smooth i ng  operators  S^,  S^, 

S^,  ...,  in  such  a way  that 

lim  = I , (5-18) 

vy.bere  I denotes  the  unit  operator;  in  other  terms,  if  Sj^f  is  the 
result  of  applying  the  smoothing  operator  to  a function  f , 

then 

lim  S^f  = f ; (5-19) 

k->“> 

SO  for  increasing  k , the  smoothing  becomes  less  and  less.  (There 
is  hardly  any  danger  to  confuse  the  smoothing  operator  with 

the  boundary  surface  S !) 

To  have  a simple  example,  let  us  take  a harmonic  function  f 
which  is  developed  into  a series  of  spherical  harmonics ; 


f 


I Y„(6,X) 


n=o  r 


n+  1 


(5-20) 


and  take 


lo+kY  (e,x) 

I " 


n = o r 


n+  1 


9 


(5-21) 
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so  that  smoothing  is  obtained  by  truncating  the  spherical  harmonic 
series  or,  what  is  the  same,  truncation  of  the  spectrum.  For 
instance , 


S f = 


10  Y^(e,A) 
n=or" 


.n+  1 


9 


11  Y (0  ,X) 

f = I " 


n=o  r 


n+  1 


> 


etc.;  it  is  clear  that  (5-19)  is  satisfied. 

The  reason  for  introducing  a smoothing  in  Mol odensky ' s 
problem  is  the  following.  It  is  a well-known  difficulty  with  many 
higher  order  solutions  that  the  higher  order  terms  are  getting 
rougher  and  rougher.  This  is  the  case  if  an  iteration  involves 
differentiation:  the  derivative  is  almost  always  less  smooth  than 
the  original  function.  A case  in  point  is  the  well-known  Molodensky 
series  for  solving  the  simple  (linear)  Molodensky  problem.  The 
calculation  of  higher  order  terms  proceeds  through  an  iteration 
which  involves  successive  differentiation,  with  the  result,  that 
higher  order  terms  can  no  longer  meaningfully  be  computed  from 
empirical  gravity  data. 

The  Hbrmander  iteration  for  the  non-linear  Molodensky 
problem  has  to  struggle  with  a similar  difficulty:  assume  we  have 
an  approximate  solution  for  the  potential  W(x;e)  according  to 
(3-34).  The  calculation  of  the  isozenithal  vector  m , which  is 
required  for  the  boundary  condition  (3-60),  requires  the  calculation 
of  the  matrix  (3-43),  involving  two  differentiations.  So  we  "lose" 
two  derivatives  at  each  iteration,  the  functions  involved  will  get 
rougher  and  rougher,  and  the  iteration  will  probably  “blow  up". 

Therefore,  we  must  counteract  this  loss  of  derivatives  by 
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a suitable  smoothing,  taking  care,  however,  to  successively  reduce 
the  amount  of  smoothing  so  that  (5-19)  is  satisfied:  otherwise  we 
would  not,  in  the  limit,  obtain  the  right  result. 

Furthermore,  the  smoothing  operator  gives  us  a natural 

way  of  obtaining  the  terms  f^^  in  (5-5);  we  put 


fo  - - C5-22) 

'k  “ - S^.jf  , (5-23) 

k = 1 , 2 , 3 , ... 

In  fact, 

- [So  * + (S,  - S,)  + + (S„  - S„  )]  f 

= SJ  , (5-24) 

and 


^ (5-25) 

o 

by  (5-19),  so  that  (5-5)  is  satisfied. 

Now  our  problem  is  to  modify  the  iteration  equations  (5-9) 
through  (5-12)  to  take  smoothing  into  account.  First  of  all,  take 
equation  (5-11).  It  is  easily  seen  to  correspond  to  the  solution 
of  the  linearized  problem  in  the  k-th  step,  as  tj;  involves  the 
linearized  operator  $ ' through  (5-13).  Thus  it  is  there  where  loss 
of  derivatives  occurs,  and  smoothing  is  appropriate.  Hence,  instead  of 
computing  at  Uj^  , we  compute  it  at 


(5-26) 


1 


L 
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and  (5-11)  is  replaced  by 


(5-27) 


Next,  consider  (5-10).  The  terra  is  already  sraoothed,  being 

defined  by  (5-23),  but  is  not,  being  given  by  (5-9)  which 

involves  $(u  ) and  ^(u  ) . The  obvious  thing  would  be  to  replace 

iC  iC  1 

e by  S e,  . , but  this  will  not  work,  as  we  shall  see. 

k-l  k.  k-l 

The  right  answer  it  to  put 


"k-i  = ^k 


’k-  1 


(5-28) 


where 


k-l 

'k  ^ ^ ^k 

j=o 


(5-29) 


represents  the  sura  of  all  errors  before  the  k-th  step,  and  to 
replace  ^ in  (5-10)  by 


^k^k  " ^k-i^k-i 


obtaining 


(5-30) 


Collecting  these  results,  we  replace  (5-9)  to  (5-12)  by 


'k-l  ■ *<“k>  - ‘'“k-l’  - '’k-l 


(5-31) 
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^ = 

(5-32) 

"k  = 

5k"k  • 

(5-33) 

«k  = 

♦(Vk)6k  ■ 

(5-34) 

"k+1 

' ''k  * ®k  • 

(5-35) 

This  is  the  Nash-Hormander  iteration  scheme  with  smoothing. 

Will  this  scheme  converge  to  the  right  solution?  A necessary 
condition  (unfortunately  it  is  by  no  means  sufficient)  is  obtained 
by  expressing  from  (5-31); 

- "("k-iI  ■ \-l  * 'k-1  ■ (5-36) 

replacing  k-1  by  k : 

and  summing  from  0 to  n : 

$(.Un^i)  - $(u^)  = . (5-38) 

O O 


By  (5-4)  this  should  tend  to  f if  n » 
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n n 

= Ifv  - S E 


■k  ‘-'k. 
o o 


n n 


(5-40) 


Now,  by  (5-29) 


n-  1 

S E = S J 

n n n k 

k=o 


which  because  of  (5-18),  can  be  expected  to  tend  to 

oo 

J'k  • 

O 

so  that,  with  (5-25),  condition  (5-39)  will  hopefully  be  satisfied. 

(It  is  easily  seen  that  (5-39)  would  not  be  satisfied  if  we 
had  put  ^k^k-i  which  was  our  first  choice  for  modifying 

(5-10)). 

Equations  (5-31)  through  (5-35)  are  the  basic  formulas  for 
our  iteration.  To  start,  we  put 

h = f 


(5-41) 


S u 

o o 


i>{\/  ) h 
o'  o 


U +6 
o o 


next,  we  put  k = 1 in  (5-31)  through  (5-35),  then  k = 2 and 
so  on . 

Hbrmander  writes  these  equations  in  a somewhat  different 
form,  introducing  a parameter  e (9^<9<®)  which  has  two  tasks: 


= X(e^)  (5-43) 


I with 

1 

i = x(0^)  , u = 1 im  x(e)  ; 

I 

i 2.  to  serve  as  a parameter  for  the  smoothing  operator, 

I such  that 

; . (5-44) 

In  both  cases,  0j^(  k = 0 , 1 , 2 , . . . ) represents  a discrete  sequence 
of  values  of  the  parameter  0 
A choice  for  0 is 


9^  and  V denoting  some  large  constants.  Then  the  difference 
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i s 


(5-47) 


and  can  thus  be  made  as  small  as  any  desired  power  of  e 
The  difference  quotient 


-1 


u,  . -u, 

k+ 1 k 


(5-48) 


can,  therefore,  be  made  as  close  as  desired  to  a derivative,  so 
that  the  discretization  error  becomes  very  small  as  k ^ “ 

Hbrmander  writes  the  basic  iteration  scheme  in  a different 
form,  which  explicitly  involves  the  parameter  e . Using  (5-23) 
and  (5-48)  and  putting 


"k  = \®k  * 


h = A g , 
k k^k  ’ 


we  transform  (5-31)  through  (5-35)  into 


g = a"  ^ 

'^k  k 


,)(f  - ,)  - A,  S^e,  : 

k k-1  ' k-l'  k-1  k k- 1 


\ ■ \“k  ■ 


“k  = ♦(\)9k  • 


U,  , = U,  + A,  U,  , 

k+1  k k k 


Here  e'  denotes  the  error  due  to  smoothing,  and  e"  expresses 

K,  -K. 


the  error  due  to  nonlinearity.  The  quantity 

E,  = (5-50) 

3 = 0 

represents  again  the  sum  of  all  errors  before  the  k-th  step. 

The  smoothing  operator  Sg  used  by  Hormander  is  rather 
similar  to  the  simple  example  (5-21).  It  also  amounts  to  a 
truncation  of  the  spectrum.  In  (5-21)  we  had  a discrete  spectrum, 
n taking  the  values  0,1,2,,..,  and  the  operator  was  also 

defined  only  for  integers  k . Now,  however,  we  have  a parameter 
9 that  runs  continuously  from  to  “ , and  Sg  has  to  be 

defined  for  continuous  e . Therefore,  one  considers  functions 
defined  on  a plane,  or  generally  in  , which  have  a continuous 

spectrum  (2  or  n-d imens i ona 1 Fourier  transform).  Now  the  smoothing 
simply  consists  in  the  application  of  a low-pass  filter,  which 
suppresses  higher  frequencies,  Teaving  lower  frequences  unchanged. 
In  other  words,  the  spectrum  of  the  function  to  be  smoothed  is 
multiplied  by  a function 

Hg(u))  = H(^)  , (5-51) 


where  H(a))  is  a symmetr i c function  of  the  frequency  w , which 
is  1 for  |a)|  < Wj  and  0 for  j w | > ; between  and  , 

H(u))  is  interpolated  by  an  infinitely  differentiable  function. 
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H(a.) 


Figure  8 

Smoothing  by  a low-pass  filter 


so  that  H(u)  is  everywhere  infinitely  often  differentiable.  This 
is  shown,  in  one  dimension,  by  Fig. 8;  the  generalization  to  higher 
dimensions  is  straightforward. 

Equivalent  to  a multiplication  of  the  spectrum  by  (5-51)  is, 
of  course,  a convolution  of  the  function  to  be  smoothed  with  a 
function  whose  Fourier  transform  is  (5-51);  cf.  ( Pa pou  1 i s , 1968 , p . 74 ) . 

The  smoothing  just  described  can  be  applied,  with  slight 
modifications,  to  functions  with  support  in  a compact  set  K , that 
• , to  functions  which  are  zero  outside  K . This  finally  gives  the 
ibility  to  extend  the  smoothing  operator  to  functions 

» npd  on  a compact  manifold,  for  instance,  on  the  earth's  surface 
»»‘nuroid.  This  is  done  by  covering  the  manifold  by  several 
' patches  and  splitting  up  the  given  function  into  a sum 

I 

Kith  compact  support  by  a "partition  of  the  unity"; 

i 
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t 

cf.  (Loomis  and  Sternberg , 1968 , p . 405 ) . 

In  view  of  smoothing  it  has  been  convenient  to  let  e run 
from  0 to  “ . This  contrasts  to  sec. 3,  where  we  had  0 < 0 < 1 . 

o = = 

It  is,  however,  easy  to  transform  a parameter  in  fo,lJ  to  one  in 
instance,  by  the  transformation 

0 = 0^  + tan(-^  Tr0  ' ) . (5-52) 


6 . Existence  and  Uniqueness  in  the 
Nonlinear  Molodensky  Problem 

The  iteration  procedure  outlined  in  the  preceding  section 
is  now  applied  to  Molodensky' s problem;  now 

u = S (6-1) 

represents  the  earth's  surface  to  be  determined,  and 

t(u)  = g (6-2) 

is  the  given  gravity  vector  on  S ; cf.  sec.l. 

The  crucial  problem  is,  of  course,  to  prove  convergence  of 
the  iteration  scheme  for  conditions  applicable  to  Molodensky's 
problem.  This  is  so  difficult  and  complicated  that  it  cannot  be 
done  here;  the  interested  reader  is  referred  to  (Hdrmander , 1975)  . 

We  shall  restrict  ourselves  to  stating  the  main  results. 
First  we  must  say  a few  words  a‘'''ut  the  norms  which  have  been  used 
in  this  problem.  A norm  is  needed  to  characterize,  in  some  sense, 
the  "size"  of  the  function  under  consideration;  it  can  be  considered 
as  a generalization  of  the  norm  of  a vector,  e.g. 
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||a||  = (a^  + (6-3) 

for  a vector  a = (a  , a , a ) 

1 2 3 

A well-known  norm  of  a continuous  function  f(x)  defined 
in  a convex  compact  set  B in  is 

llfll^  = suplf(x)l  , (6-4) 

xeB 

where  "sup"  denotes  the  supremum  (which,  for  continuous  functions, 
equals  the  maximum  value).  The  functions  f for  which  the  norm 
(6-4)  is  finite  form  a Banach  space;  cf.  (Loomis  and  Sternberg, 
1968, p. 218).  This  space  will  be  denoted  here  by  H° 

An  appropriate  norm  for  functions  that  are  differentiable 
as  well  as  continuous  is 


llfll  = sup  jf  (x)|  + sup 

xeB  x6B 


3X. 


(6-5) 


where  9f/3x,  denotes  any  partial  derivative.  The  functions  with 
fintie  norm  (6-5)  form  a Banach  space 

Now  it  is  of  basic  importance  for  Molodensky's  problem,  as 
for  many  problems  in  potential  theory,  to  define  norms  IKII^  ^or 
0<a<l  , that  is,  spaces  H“  intermediate  between  H°  and 

For  this  purpose  we  consider  continuous  functions  that  satisfy  a 
Holder  condition  with  exponent  a ; they  are  functions  for  which 


sup 

X ,yeB 


is  finite.  These  functions  form  a Banach  space  H“  ; the  norm  is 
given  by 
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llfll  = suplf(x)j  + sup  I I , (6-7) 

x€-B  X , yeB  I x-y  I 

It  can  be  shown  that 

H°  D 3 (6-8) 

for 

0<a^<a2<l;  {6"9) 

that  is,  there  are  more  functions  in  H°  than  in  , more 

functions  in  than  in  8“^  , and  more  functions  in  than 

in  . So  satisfying  a Holder  condition  with  exponent  a is  a 

stronger  condition  than  mere  continuity  and  weaker  than  differentia- 
bility. 

We  may  also  consider  a Holder  condition  (6-6)  with  exponent 
a = 1 ; it  is  seen  that  this  is  almost  (although  not  completely) 

the  same  os  differentiability.  In  fact,  we  shall  use  H^  for  the 
space  of  functions  satisfying  a Holder  condition  with  a = 1 rather 
than  for  functions  with  finite  norm  (6-5). 

So  far,  we  have  defined  spaces  H^  for  0<a<l  . For  a 

>1  we  proceed  as  follows.  Let  k be  a positive  integer  such  that 
k<a<k+l  (for  instance,  for  a = 5.75  there  is  k = 5 ).  Denote 
by  D f any  derivative  of  k-th  order  (for  instance. 


3 

3 Xj3  X2 

for  k = 5 ).  Then  the  norm  ||f||^  is  defined  as 


L 
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\\f\\^  = suplf(x)[  + sup  I.p''.nx)-D|^^(y)[  . (6-10) 

x«B  x,y€-B  I x-y| 

It  is  clear  that  (6-7)  is  a special  case  of  (6-10)  with  k = 0 
In  other  words,  the  space  consists  of  continuous 

functions  which  are  k times  differentiable  and  whose  k-th 
derivatives  satisfy  a Holder  condition  with  exponent  a-k<l 

Reformulation  of  the  Molodensky  Problem.-  Hormander  has 
reformulated  Molodensky's  problem  in  the  following  way:  to  determine 
a closed  surface  S in  , which  is  a one-to-one  image  of  the 

unit  sphere,  from  given  values  "g  and  W , such  that  the  following 
conditions  are  satisfied: 


3 


W = WoS  + ya  A 

1 ^ 3 

> 

(6-11) 

g = gcS  = W'oS 

> 

(6-12) 

W(x)  = V(x)  + ^ 

<u"(x2  + x^)  . 

(6-13) 

aV  = 0 outside 

s . 

(6-14) 

V(X)  = 1 + 

o(-^)  . 

1 xl  ^ 

(6-15) 

Uniqueness  of  the  solution  is  achieved  by  postulating  that  the 
harmon i c function  V(x)  , which  represents  the  external  gravitational 
potential,  contains  no  first-degree  spherical  harmonics;  this  is 
expressed  by  (6- 15 ) , 

jxl  = r (6-16) 

denoting  the  radius  vector.  The  second  term  on  the  right-hand  side 
of  (6-13)  denotes,  of  course,  the  centrifugal  force  potential,  w 
being  the  angular  velocity  of  the  earth's  rotation. 
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The  new  feature  is  (6-11)  instead  of  simply  taking  ¥ = WoS  , 
which  would  be  the  restriction  of  W to  S .In  the  modified 
expression,  the  a^  are  three  constants  to  be  determined,  and  the 
are  three  suitably  assumed  functions.  The  purpose  of  adding  the 
linear  combination  ensure  the  solvability  of  Molodensky's 

problem  for  arbitrar^'  boundary  data. 

This  is  to  be  understood  as  follows.  Assume  the  earth's  sur- 
face S to  be  known,  and  consider  the  given  function  ¥ = WoS 
We  can  now  solve  the  exterior  Dirichlet  problem 

= 0 , V ->  «>  for  r ->■  oo  , 

VoS  = ¥ - a)2(s2  ^ 52)  ^ (6-17) 

which  gives,  for  every  data  function  ¥ , a unique  solution  V(x) 

For  arbitrary  data  ¥ , the  spatial  function  V(x)  will,  in 

general,  contain  spherical  harmonics  of  first  degree,  contrary  to  the 
condition  (6-15)  . 

Let  now  (6-11)  be  used  instead  of  ¥ = WoS  . Then  (6*17) 
is  replaced  by 

VoS  = ¥ - -^  a)^(sj  + S^)  ' la.A^  . (6-18) 

Again,  the  solution  of  the  exterior  problem  with  boundary  data  Vo S 
gives  a unique  solution  V(x)  which  will  contain  three  linearly 
independent  spherical  harmoni cs  of  the  first  degree.  Now,  however, 
the  three  constants  a^  can  be  chosen  in  such  a way  that  these 
three  first-degree  harmonics  vanish. 

This  is  readily  seen  to  lead  to  three  linear  equations  for 
the  three  unknowns  a^,  a^,  a^  . These  equations  will  have  a unique 
solution,  provided  the  boundary-value  problem 
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3 

AV  = 0 . VoS  = p.A.  (6-19) 

has  a solution  v(x)  which  contains  three  linearly  independent 
forst-degree  harmonics 

(cosine  cosx  + cosine  sinx  + C2Cose)/r^  . (6-20) 


In  fact,  c^  can  then  be  chosen  to  be  equal  to  the  c^  of  the 
boundary-value  problem  (6-17),  and  similar  for  c^  and  c^  . Since 
the  problem  (6-18)  is  the  difference  of  problems  (6-17)  and  (6-19), 
all  first-degree  harmonics  will  cancel  in  the  solution  of  (6-18). 

To  achieve  this,  we  may  select  the  functions  as  follows: 


SI  ne  cosx 


)oI 


s 1 ne  s 1 nx 


)oI 


/ cose . V 
( — 


(6-21) 


If  the  earth's  surface  S coincided  with  the  telluroid  I , then 
the  spatial  functions  corresponding  to  the  A^  by  a solution  of  the 
exterior  Dirichlet  problem  would  be  the  first-degree  harmonics 


sine cosx 


, etc . 


(6-22) 


themselves,  for  which  the  condition  (6-20)  is  certainly  satisfied 
(with  c^  = a ) . If  S does  not  deviate  too  much  from  T , then 
(6-20)  is  still  satisfied  b e c a u 

A^($,A)  are  now  considered  as  functions  on  S instead  of  I (now, 
of  course,  c^  / a^  in  general).  More  generally,  if  the  A^(i>,a) 
as  given  by  (6-21)  independently  of  the  parameter  e , are  considered 

as  functions  of  e with  S(0^)  = I (telluroid)  and  S(oo)  = s 
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(earth's  surface),  then 

a.  = a.(e)  (6-23) 

will  be  functions  of  0 

The  linearization  of  the  modified  Molodensky  problem  defined 
by  (6-11)  through  (6-15)2is  now  done  as  in  sec. 3.  The  only  difference 
is  the  additional  term 
(3-60)  is  replaced  by 

• • • • J 

(V  + <m,grad  V)o)]  = IT  + <mo)],g'>  - ^a.A,  , (6-24) 

1 ^ ^ 


so  that  the  boundary  condition 


where  a.  is  the  derivative  of  (6-23)  with  respect  to  9 . Now  for 

arbitrary  boundary  data  w and  g there  exist  unique  constants 
a,,  a * a , so  that  the  corresponding  boundary-val ue  problem  with 
= 0 outside  )|(e)  admits  a unique  solution  V(x)  that  does  not 
contain  spherical  harmonics  of  first  degree. 

To  the  Molodensky  problem  reformulated  in  this  way,  the 
Nash-Hbrmander  method  is  applied.  Convergence  of  this  method  can 
be  proved  using  various  Holder  norms  and  estimates  for  the  linear 
problem  and  also  for  the  nonlinearity.  All  this  is  extremely  difficult 
and  laborious.  Finally  one  obtains  the 

Theorem  of  Hbrmander.  Let  e > 0 , then; 

(1)  Fop  all  W and  g in  a H neighborhood  of  W 

and  g^  j then  the  modified  Molodensky  problem  defined  by  equations 
(6-11)  to  (6-15)  has  a solution  S close  to  S = in  and 

■L__  ,.2—  ° ^ 

(2)  If  W and  g are  in  for  some  a>2+E  which  is  not 

an  integer,  then  SeU^  . of 

(3)  One  can  find  a neighborhood  which  cannot  contain 


two  solutions  of  the  problem. 
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Let  us  look  at  this  theorem  more  closely  and  explain  its 
meaning,  A neighborhood  of  = ¥(u;e^)  consists  of  all 

functions  ¥ for  which 

. (6-25) 

where  6 is  sufficiently  small  and  the  norm  is  defined  by  (6-10)  with 
a = 2 + e . Smallness  of  this  norm  implies  that  not  only  the  maximum 

deviation  of  ¥ from  ¥ , 

o 

max j¥  - ¥ I 
'o' 

issmalljbutalsothat  ^ 

max|¥'  - ¥^ I 

and 

maxi¥"  - ¥" | 

is  small,  so  that  not  only  ¥ must  be  close  to  ¥^  , but  also  the 

first  and  the  second  derivatives  of  ¥ must  be  close  to  those  of 
¥ . In  addition  to  this,  something  more  is  required.  If  e = 1 , 

O 

then  also  closeness  of  the  third  derivatives  must  hold;  if  0<e<l  , 

the  intermediate  Holder  condition  is  stronger  than  mere  closeness  of 
the  second  and  weaker  than  closeness  ot  the  third  derivatives:  the 
difference  of  the  second  derivatives  must  satisfy  a Holder  condition. 

Closeness  of  the  telluroid  I and  the  earth's  surface  S 
in  H^'*’^  means  that  the  maximum  deviation  of  the  surface  is  small 
and  that,  in  addition,  slopes  (first  derivatives)  and  curvatures 
(second  derivatives)  are  also  similar  for  S and  I ; in  addition, 
there  is  a Holder  condition  for  the  difference  of  the  second  deriva- 
tives. 


I 
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Part  (1)  of  Hormander's  theorem  asserts  the  exi stence  of  a 
solution  provided  we  have  a good  approximation  ^ for  the  earth's 
surface  S not  only  with  respect  to  the  maximum  deviation  between 
S and  ^ , but  also  with  respect  to  first  and  second  derivatives 

(plus  a Holder  condition),  and  also  a good  approximation  to  potential 
and  gravity. 

This  condition  is  obviously  very  strong.  If  one  uses  an 
ellipsoidal  reference  field  and  the  telluroid  according  to  the  usual 
definition,  then  the  actual  gravity  field  and  the  earth's  surface 
probably  fall  short  of  this  condition.  It  is,  however,  not  required 
that  the  initial  approximations  for  S and  W satisfy  this  condition; 
it  would  be  sufficient  if  any  intermediate  approximation  would  meet 
it  (because  then  this  intermediate  approximation  could  be  considered 
as  the  starting  point).  Still  it  is  not  clear  even  then  whether 
Hormander's  theorem  could  be  applied  to  the  actual  earth. 

Part  (2)  of  the  theorem  assures  that  the  resulting  surface 
S will  be  as  smooth  as  the  data:  if  the  data  are  n times  differen- 
tiable and  if  the  n-th  derivatives  satisfy  a Holder  condition,  then 
the  same  will  hold  true  for  S 


Part  (3)  ensures  uniqueness  but  under  an  even  stronger  con- 
dition,  (H  neighborhood)  than  for  the  existence  theorem  of 
Part  (1)  (H^^^  neighborhood).  However,  Hbrmander  thinks  it  highly 
probable  that  H^^^  could  be  replaced  by  H^^^  , so  that  the 

condition  for  uniqueness  would  be  the  same  as  for  existence. 

In  Part  2,  integer  values  of  e are  excluded;  this  reflects 
the  well-known  fact  that  Holder  conditions  with  e ^ 0 are  essential 
in  potenti al -theoretical  considerations.  In  Parts  1 and  3,  also 
integer  e are  admitted. 

In  conclusion  we  may  say  that  Hormander's  theorem,  although 
not  directly  applicable  to  the  real  earth,  gives  the  first  mathe- 


matically exact  results  on  existence  and  uniqueness  for  Molodensky's 
problem  and  is  thus  of  fundamental  importance. 
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7 . The  Gravity  Space  Approach 


Recently,  F.  Sansd  (1977  a,b,c,d)  has  given  a completely 
different  approach  to  the  nonlinear  Molodensky  problem.  The  idea  is 
to  use  the  three  cartesian  components  of  actual  gravity,  g^ , g^ , 93 
as  new  curvilinear  coordinates,  instead  of  the  cartesian  coordinates 


X , X , 
1 > 2 ’ 

of  the 


themselves.  Thus  the  potential 


W = W(£)  = W(g^  , g^  , g^ 


becomes  a function 


(7-1) 


On  the  physical  earth's  surface  S , the  three  components  g^  of 
the  vector  £ are  given,  as  well  as  the  potential  W ; therefore 
the  three  curvilinear  coordinates  g^,  g^ , 93  of  each  point  of  the 
surface  S are  known,  or  S is  a known  surface  jf  expressed  in 
terms  of  coordinates  g^ 

There  are  two  ways  of  interpreting  g : either  they  may  be 

i 

considered  as  curvilinear  coordinates  in  ordinary  space,  or  as 
cartesian  coordinates  in  an  auxiliary  space,  called  gravity  space. 
Using  the  second  interpretation,  we  may  say  that  S becomes  a known 
surface  in  gravity  space  or  the  free  boundary-value  problem  is 
transformed  into  a fixed  boundary-val ue  problem.  The  simplification 
which  is  achieved  in  this  way  for  the  nonlinear  Molodensky  problem 
is  so  decisive  that  inconveniences  and  difficulties  arising  with 
this  indirect  approach  are  more  than  compensated  as  far  as  theoreti- 
cal investigations  on  existence  and  uniqueness  of  the  solution  are 
concerned . 


The  main  inconvenience  is  that  a one-to-one  correspondence 
between  cartesian  coordinates  x.  and  g.  outside  and  on  S , 
which  is  the  region  in  which  we  work,  is  possible  only  if  the  earth 
is  nonrotating.  To  see  this,  consider  ellipsoidal  gravity  along 
a radius  vector  in  the  equatorial  plane.  As  the  height  increases. 


1 
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Y first  decreases  but  then  it  increases  again  because  the  centri- 
fugal force  becomes  dominant.  So  at  a certain  elevation,  y will  be 
the  same  as  on  the  ground,  which  violates  a one-to-one  correspondence 
between  gravity  vector  and  position.  For  the  actual  gravity  field 
the  situation  is  similar  as  in  the  ellipsoidal  case. 

If  the  earth  is  considered  as  nonrotating,  then  the  corre- 
spondence between  gravity  and  position  is  seen  to  be  one-to-one 
(provided  the  Marussi  condition  holds,  see  below).  In  other  terms, 
the  correspondence  is  unique  if  we  work  with  the  gravitational 
potential  V and  the  gravitational  vector  grad  V instead  of  the 
gravity  potential  W and  the  gravity  vector  g = grad  W 

It  is,  of  course,  clear  that  only  W and  grad  W (including 
the  centrifugal  force)  are  directly  measurable.  However,  the  effect 
of  centrifugal  force  can  be  calculated  with  sufficient  accuracy  on 
the  basis  of  our  current  knowledge  of  the  earth's  surface  (the  error 
in  the  centrifugal  force  is  less  than  + 0.005  mgal  for  a position 
error  of  + 10  meters),  and  subtracted  from  W and  grad  W to  give 
their  gravitational  counterparts  V and  grad  V . Therefore,  Sanso's 
boundary-value  problem,  which  uses  gravitation  instead  of  gravity, 
is  practically  as  meaningful  as  the  original  Molodensky  problem. 

In  the  sequel  we  shall  thus  work  with  the  gravitational  po- 
tential V , which  is  a harmonic  function,  and  take  ^ as 

£ = grad  V , (7-2) 

so  that  £ is  the  vector  of  gravitation  rather  than  gravity.  We 
shall,  however,  continue  to  call  g , even  if  defined  by  (7-2),  the 
gravity  vector,  to  be  in  agreement  with  Sansb's  terminology  and 
with  the  term  "gravity  space"  (this  is  consistent  with  current 
terminology  if  we  consider  the  earth  as  nonrotating). 


— J 
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We  can  thus  reformulate  Mol odensky ' s problem  in  terms  of 
V as  follows;  to  find  a function  V ( x)  which  is  harmonic  outside 
a unknown  closed  surface  S , 


aV  = 0 


(7-3) 


and  which,  together  with  its  gradient,  assumes  on  S the  given 
boundary  values 


VoS  = V’{u) 

(grad  V)oS  = £(u) 


(7-4) 

(7-5) 


where 


u = («>,a) 


(7-6) 


as  in  sec  . 3 . 


We  now  introduce  the  components  g^  of  g as  new  coordi- 
nates, which  are  functions  of  the  rectangular  coordinates  (x^,  x^, 


X3)  : 


9i  = • 


(7-7) 


If  his  transformation  is  to  have  an  inverse. 


Xj  = x,(3^)  . 


(7-8) 


then  the  Jacobian  determinant 


det 
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must  be  nonzero  everywhere  on  and  outside  S . Since 


n - 

" 3X . ’ 

1 


this  condition  is 


(7-9) 


det 


3^V 


3 X . 3 X . 
1 3 


i 0 


(7-10) 


which  is  nothing  else  than  the  Marussi  condition;  this  condition 
(restricted  to  S ) has  al ready  been  used  before;  cf.eq.  (3-48).  It 
will  be  assumed  that  the  Marussi  condition  is  satisfied  everywhere 
outside  and  on  S 

Now  the  potential  V becomes  a function  of  the  vector  £ : 

V = V(£)  = V(gj,  g^.  g3)  (7-11) 

As  we  have  mentioned,  this  would  reduce  Molodensky's  problem  to  a 
fixed  boundary-value  problem  (actually  a Dirichlet  problem)  in 
gravity  space.  Since  V as  a function  of  £ satisfies  a linear 
partial  differential  equation  of  second  order,  which  is  Laplace's 
equation  aV  = 0 , it  does  the  same  as  a function  of  £ since 

the  transformation  (7-7)  or  (7-8)  transforms  Laplace  equation  into 
another  linear  second-order  partial  differential  equation.  However, 
since  the  transformation  (7-8)  is  actually  unknown,  the  coefficients 
of  this  differential  equations  are  noc  known,  and  therefore  this 
approach  appears  hopeless. 

Sansb  has  found  an  ingenious  way  out  of  this  difficulty  by 
transforming  not  only  the  coordinates  but  also  the  potential, 
introducing  an  adjoint  potential 
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y = x.£  - V = X.  g - V ; 


(7-12) 


this  is  a Legendre  transformation  familiar  from  other  fields  (ordi 
nary  differential  equations,  analytical  mechanics,  thermodynamics, 
etc.) . 

Differentiating  (7-12)  with  respect  to  g^  we  get 


g„  + X,- 


av 


3g.  3g.  =’k  i 3X,  9g. 

1 k 


= X 


in  view  of  (7-9)  . Thus 


ay 


(7-13) 


or 


X = gradin'  , 


(7-14) 


which  shows  a striking  symmetry  betiveen  x^  and  V on  the  one 
hand  and  g^  and  \\)  on  the  other  hand. 

Also  (7-12)  is  completely  symmetric 


V + f = x^gj^ 


(7-15) 


and  permits  to  express  one  potential  in  terms  of  the  other; 


H'  = X 


3V 
k 3X, 


V 


(7-16) 


V = g 


39^ 


- y . 


The  matrix  of  second  gradients  of  'V  , 


-'i'  |_39i39j 


(7-17) 


(7-18) 


(by  (7-13))  is  inverse  to  the  matrix  of  second  gradients  of  V : 


3x.  9X_, 
_ 1 j 


3x^  ’ 


cf.  (2-35)  and  (2-36);  that  is. 


(7-19) 


My  = 


(7-20) 


Now  Laplace's  operator 


,v  = iil . 

axf  dxt  ax, 

1 2 3 


(7-21) 


is  nothing  else  than  the  trace  Tr  of  the  matrix  , and 

Laplace's  equation  may  be  written 


Tr  M = 0 

-V 


(7-22) 


This  gives  us  a possibility  to  find  the  corresponding  partial 
differential  equation  for  ^(gj,  g^)  : by  combining  (7-20) 
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and  {1 -11)  we  get 


Tr(M;^)  = 0 


(7-23) 


On  introducing 


f . 


2 

3 4' 


ij  ag^ag^ 


(7-24) 


the  matrix 

M 

1^  becomes 

y 

11 

y 

12 

y 

13 

My 

= 

y 

21 

y 

22 

y 

23 

i 

1 

y 

L 31 

y 

32 

y 

33 

(7-25) 


Inverting  this  matrix  and  taking  the  trace  gives 


2 2 
mm  _m'‘+mu(  -m 

^11  22  12  22  33  23  11  33 


y = 0 
13 


(7-26) 


which  is  a partial  differential  equation  for  'i’(g^,  g^,  g^)  with 
known  coefficients  (all  + 1),  but  unfortunately  a nonl i near  one. 

(In  two  dimensions,  we  would  again  get  Laplace's  equation  for 

'^(91*  93)  •) 

The  basic  differential  equation  (7-26)  may  also  be  written 
in  the  form 


(Tr  M^)2  - Tr(M2)  = 0 . 


(7-27) 


which  is  verified  by  direct  calculation. 

In  gravity  space,  the  vector  £ is  the  position  vector, 
the  components  serve  as  rectangular  coordinates  and  gravity 

g serves  as  radius  vector.  In  fact,  we  have  from  (1-7) 

g^  = g cos$  cosA  , 

g^  = g cos$  sinA  , (7-28) 

93  = g sin$  , 

where  $ and  a are  the  astronomical  coordinates  (for  a non- 
rotating earth  or  after  removal  of  centrifugal  effects).  This  shows 
that  g,  <i>,  A are  nothing  else  than  spherical  polar  coordinates 
in  gravity  space.  The  derivative  3/3g  is  thus  a radial  derivative 
in  gravity  space;  we  have 

11  = 9T 

3g  3gj^  3g  ' 39^  9 
using  (7 -28)  . Thus 


11  - n 

39  ■ ^k  3gj^ 


and  (7-17)  may  be  written  as 


V = 


(7-29) 


(7-30) 


The  boundary  condition  in  gravity  space  thus  becomes 
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(9  If  ■ » (7-31) 

where  the  known  function  Y^u)  is  given  as  a function  of  the 
parameter  (7-6)  which  in  gravity  space  denotes  the  two  angular 
spherical  coordinates;  is  the  image  of  the  earth's  surface  in 

gravity  space. 

For  large  values  of  the  spatial  radius  vector 
r = = lx|  (7-32) 

we  have 

V = ^ + O(^)  , (7-33) 

r 

9 = -^  + O(l^)  , (7-34) 

r ^ r 

where 

n = GM  (7-35) 

denotes  the  product  of  the  gravitational  constant  G and  the 
earth's  mass  M ; we  have  taken  the  coordinate  origin  at  the 
earth's  center  of  mass. 

For  r -►  ® we  have  g ->■  0 , so  that  the  spatial  infinity 

corresponds  to  the  origin  in  gravity  space.  Solving  (7-34)  for  1/r  , 

-1  i 

i = y'2g2  ^ ^ (7.3g) 
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and  substituting  this  into  (7-33)  we  get 
1 i 1 

V = >i^g^  + 0(g2)  , (7-37) 

which  expresses  the  behavior  of  V as  g ^ 0 . Finally, 

i 1 1 

= - 2p2g2  + 0(g2)  . (7-38) 

which  is  verified  by  substitution  into  (7-30),  taking  (7-37)  into 
account . 

We  thus  arrive  at  the  following  formulation  of  the  geodetic 
boundary-value  problem  in  gravity  space:  to  find  the  solution  of 
the  partial  differential  equation  (7-26)  outside  with  the 

boundary  condition  (7-31)  on  ; the  earth's  surface  S will 

then  be  given  by  (7-14): 

2<oS  = (grad^'C  )oSg  , (7-39) 

where  xoS  denotes  the  position  vector  x restricted  to  the  sur- 
face S , that  is,  the  position  vector  of  any  surface  point, 
x(^,A)  . 

Since  the  direction  3/3g  is  the  direction  of  the  radius 
vector  in  gravity  space,  in  general  different  from  the  normal  to 
, we  have  an  oblique-derivative  problem  with  a known  surface 
Sg  and  a linear  boundary  condition  (7-31),  but  for  a nonlinear 
partial  differential  equation  (7-26). 
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8 . Linearization;  Comparison  to  the 
Standard  Approach 

The  linearized  equation  (2-43)  shows  a striking  formal 
analogy  with  (7-17).  To  take  a closer  look  at  this  analogy,  we  shall 
also  linearize  (7-17)  and  other  relations  in  gravity  space. 

We  shall  use  the  concept  of  the  gravimetric  telluroid 
explained  in  sec. 2;  there  is  a one-to-one  correspondence  between 

I 

the  points  P of  the  earth's  surface  S and  Q of  the  gravimetric  j 

telluroid  S by  postulating  j 

■ 

y,(Q)  = g.(P)  . (8-1) 


that  is,  the  normal  gravity  vector  at  Q is  to  be  equal  to  the 
actual  gravity  vector  at  P ; cf.  Fig. 2. 

As  always  in  the  gravity  space  approach,  we  assume  the 
earth  as  nonrotating  or,  which  is  the  same,  the  potential  is  the 
gravitational  potential  V . The  normal  gravitational  potential 
will  be  denoted  by  V . Then  the  disturbing  potential  T is 

T = V - V ; (8-2) 

it  is  the  same  as  in  the  usual  definition  T = W - U since  the 
centrifugal  potential  cancels  in  the  difference. 


The  adjoint  potentials  corresponding  to  V and  V are 
given  by  (7-12): 


98 


Here  we  have  been  careful  in  specifying  the  arguments.  The  gravity 
space  for  normal  gravity  is  identified  with  the  gravity  space  of 
actual  gravity:  equal  numerical  values  of  g.  at:d  y.  correspond 
to  the  same  point  in  gravity  space.  It  is,  therefore,  possible  to 
denote  the  independent  variable  in  gravity  space  simply  by  g^  , 
also  when  the  normal  potential  is  under  consideration,  for  instance, 
in  (8-4). 

Equations  (7-7)  and  (7-8)  give  the  transformation  between 
ordinary  space  and  gravity  space  for  actual  gravity.  The  correspon- 
ding transformations,  between  ordinary  space  and  gravity  space,  for 
normal  gravity  are  given  by 


9i  = Ti(Xj)  , 

(8-5) 

* 51(93)  • 

(8-6) 

In  (8-5),  g^  denote  the  coordinates  in  gravity  space,  and  T^(Xj) 
are  the  functions  which  express  normal  gravity  in  terms  of  the 
coordinates  x.  ; the  £.(g.)  in  (8-6)  are  the  inverse  functions 

j 1 j 

of  Y^(x.)  • This  will  explain  the  notation  used  in  (8-3)  and  (8-4). 

It  is  clear  now  that 

x.(g.)  = x.(P)  (8-7) 

are  the  coordinates  of  the  point  P and 

Ci(gj)  = x.(Q)  (8-8) 

are  the  coordinates  of  the  point  Q , in  view  of  (8-1);  for  the 
same  reason,  P and  Q are  mapped  into  the  same  point  in  gravity 
space : 
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(8-9) 


to  which  both  't'(g^)  and  v(g^)  in  (8-3)  and  (8-4)  refer 
Let  us  now  calculate  the  difference 


T = ¥ - '!' 


(8-10) 


which  is  the  gravity  space  equivalent  of  the  disturbing  potential 
T as  given  by  (8-2).  Subtracting  (8-3)  and  (8-4)  we  get 


- » [><3(91)]  * V [EjOj)] 

In  agreement  with  (2-15)  we  put 


(8-11) 


(8-12) 


(we  omit  the  argument  g^  : denotes  the  coordinates  of  P and 

Cj  those  of  Q ),  Now  to  V [^j(9j^)]  we  apply  Taylor's  theorem: 

V(Xj)  = V(c,  + c.) 

3V 


= ^ ft  S. 


(8-13) 


The  substitution  of  (8-12)  and  (8-13)  into  (8-11)  gives 
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T(gi)  = - V(c.)  - + V(c.) 

= - V(C,)  + V(5j) 

T{g.)  = - T [Cj(g.)]  . (8-14) 


In  geometrical  terms,  x at 
Q . 


equals  the  negative  of  T at 


We  have  thus  obtained  the  result  that  the  adjoint  potential 
of  T is  simply  the  negative  of  T . This  is  certainly  surprising 
at  first  sight,  and  it  indicates  a deep  relation  between  gravity 
space  and  ordinary  space:  gravity  space  is  not  just  an  artifice 
introduced  ad  hoc  , but  a natural  expression  of  the  mathematical 
structure  of  the  geodetic  boundary-va 1 ue  problem. 

This  will  even  become  more  evident  if  we  consider  the 
boundary  condition.  In  view  of  (8-1)  we  have 


S = 
g 


(8-15) 


so  that  the  earth's  surface  S and  the  telluroid  I are  mapped 
into  the  same  surface  S in  gravity  space. 

By  (7-29),  the  boundary  condition  (7-31)  becomes 

(g^  - ’jos,  • T(u)  . (8-16) 

The  corresponding  condition  for  the  norma  1 potential  V at  the 
telluroid  I is 


101 


at 


(8-17) 


The  subtraction  of  these  two  equations,  which  are  linear  in  'V  and 
? , gives  by  (8-10) : 


(g 


9 T 

k 3g, 


- t)oS  = T(u)  - V(u) 


(8-18) 


Now , 


V’(u)  - V(u)  = ¥(u)  - II(u) 


= W - U 
P Q 


because  TJ  refers  to  S and  H to  and  because  the  difference 
of  the  centrifugal  potentials  at  P and  at  Q is  negligibly  small, 
By  ( 2 - 3 ) this  is 


7(u)  - V(u)  = AW  . 


(8-19) 


Furthermore  t on  equals  - T on  J . Thus  (8-18)  becomes 


- g 


iL 

k ag. 


+ T = aW  , 


(8-20) 


which  is  now  a boundary  condition  on  the  telluroid  1 . The 

replacement  of  g^^  by  changes  (8-20)  only  be  second-order 

quantities,  which  are  to  be  neglected.  Thus  the  boundary  condition 
on  I finally  takes  the  form 


iL 

'k  ay, 


aW  . 


T 


(8-21) 


1 


I 
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This  is  nothing  else  than  (2-43)  with  = 0 for  the  gravi- 
metric telluroid,  and  with  (2-44).  We  thus  have  recovered  the 
fundamental  boundary  condition  of  sec. 2 via  gravity  space. 

What  about  the  differential  equation  which  t must  satisfy? 

We  could  derive  it  from  (7-26),  but  there  is  a much  simpler  way, 
using  (8-14).  In  this  equation  we  substitute  (8-5)  and  (8-6), 
obtaining 

T [Yi(x,)]  = - T(x.)  (8-22) 

Since  T satisfies  Laplace's  equation 

at  = 0 , (8-23) 

T = - T will  also  satisfy  it; 

At  = 0 -,  (8-24) 


if  T is  considered  as  a function  of  » then  the  Laplacian  is 

to  be  expressed  in  terms  of  , which  here  are  to  be  regarded  as 

curvilinear  coordinates  in  ordinary  space  related  to  the  x^  by 
(8-5).  It  is  not  difficult  to  transform  the  Laplacean  to  curvilinear 
coordinates;  cf.  ( Hot i ne , 1969 , p . 19 ) ; the  important  thing  to  note  is 
that  it  does  not  have  the  "cartesian"  form: 


(8-25) 


Thus,  as  far  as  the  linear  problem  goes,  the  gravity  space 
approach  simply  amounts  to  the  use  of  curvilinear  coordinates 
in  ordinary  space.  It  is,  therefore,  not  essentially  different  from 
the  usual  approach  outlined  in  sec. 2;  it  is  even  less  general  as  it 
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supposes  a nonrotating  earth.  The  situation  is  quite  different  for 
the  noni i near  problem  where  the  gravity  space  approach  introduces 
essentially  new  features  and  a considerable  simplification. 

Different  as  the  ordinary  approach  and  the  use  of  gravity 
space  are,  the  linearized  problem  is  the  same  in  both  methods.  This 
is  practically  important  because  the  linearized  Molodensky  problem 
is  probably  sufficient  for  all  present  applications,  as  we  have 
pointed  out  at  the  end  of  sec. 2. 

Even  for  the  linear  problem,  however,  the  gravity  space 
approach  provides  a deeper  insight  into  the  problem;  in  particular, 
the  structure  of  the  operator  that  acts  on  T in  (2-43), 


T 


aT 

''i  377 


(8-26) 


is  interpreted  by  the  relation  between  potential  and  adjoint 
potential  as  expressed  by  (7-17). 

Spherical  Approximation.-  Let  us  finally  introduce  a 
spherically  symmetric  normal  potential;  this  corresponds  to  the 
"spherical  approximation"  outlined  in  sec. 2. 

For  a spherically  symmetric  mass  configuration  we  have 


cf.  (7-33).  By  a simple  change  of  scale  of  length  and  without  loss 
of  generality  we  can  make  y = 1 , obtaining 


V 


_1 

r 


(8-27) 
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Y = 


(8-28) 


with 


(8-29) 


It  is,  therefore,  possible  to  express  the  in  terms  of  by 

2 

^ (8-30) 


X.  = - Y Yi 


in  the  case  of  a spherically  symmetric  mass  configuration,  cartesian 
coordinates  x^  and  gravimetric  coordinates  are  thus  related 

in  a simple  way. 

Another  possibility  to  convert  gravimetric  coordinates 
into  cartesian  coordinates,  denoted  by  y^  , is  by  putting 


Yi  = Y 


(8-31) 


These  coordinates  y^  can  be  interpreted  in  the  following  way.  Let 
us  consider  an  inversion  in  the  unit  sphere  r = 1 ; cf.  (Kellogg, 

1929, p. 231).  This  inversion  transforms  a point  with  coordinates 
x^  into  a point  with  coordinates  xj^  given  by 
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x:  = ^ X.  (8-32) 

r 

the  inverse  transformation  being 

X.  = x!  with  r'^  = x>’  . (8-33) 

1 ^ , 2 1 k k ' ' 

On  substituting  (8-30)  and  comparing  the  result  with  (8-31)  we  see 
that 


= - x;  , (8-34) 

so  that,  apart  from  the  sign,  are  the  cartesian  coordinates  of 

the  image  of  the  point  x^  under  an  inversion  in  the  unit  sphere. 

The  corresponding  transformation  of  harmonic  functions  is 
called  a Kelvin  transformation  ( i b i d . , p . 2 3 2 ) . The  basic  principle 
is  that  if  U(x^)  is  a harmonic  function  of  x^  in  a domain  T , 
then 


X ' 

V(x:)  = ^ U(^)  (8-35) 

is  a harmonic  function  of  x!  in  the  domain  T'  into  which  T is 

1 

carried  by  the  inversion. 

So  far,  we  have  interpreted  this  transformation  as  a 
point  transformation,  which  transforms  a point  P(x^)  into  a point 
P'(x^)  , the  coordinates  x^  and  x^'  referring  to  the  same 

cartesian  coordinate  system.  We  may , however,  interpret  it  also  as 
a coordinate  transformation,  by  which  the  same  point  in  space  i '• 
referred  to  different  coordinate  systems  x^  and  x!^  . Then  the 

Kelvin  transformation  implies  that  if  U(x^)  satisfies  Laplace's 
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equation  in  "cartesian"  form  using  x. 


1 u = ^ = 0 


3X‘ 


3X' 


3x! 


(8-36) 


then  the  function  (8-35)  satisfies  Laplace's  equations  in  cartesian 
form  using  x : 


A V= 


3X 


3X 


. = 0 
3X^^ 


(8-37) 


In  view  of  (8-34),  Laplace's  operator  will  then  have  cartesian  form 

also  in  coordinates  y : 

i 


A V = 0 
y 


(8-38) 


The  symbol  a , A , etc.  will  be  reserved  for  Laplace's  operator 
X y 

in  cartesian  f orm . 

Let  us  now  apply  these  considerations  to  the  present  problem. 

We  have  seen  that  T(x^)  satisfies  Laplace's  equation  a^T  = 0 ; 

cf.  (8-23).  t(y^)  also  satisfies  Laplace's  equation  (8-24),  but 

not  in  cartesian  form  (8-25).  If  in  t(y  ) we  introduce  new  coordinates 

i 

y^  , defined  by  (8-31),  putting 


= y^y  with  y^  = yj^yj^  . 


(8-39) 


then  the  new  function 


♦(y,)  = y 


T 

y 


(8-40) 


107 


will  satisfy 

= 0 (8-41) 

because  of  (8-35)  with  U = - T , V = <()  and  y.  = - x!  , since  T 
satisfies  aT=0 

X 

Also  the  function  v defined  by 

v(yi)  = - 4>)  (8-42) 

is  harmonic: 

A V = 0 . (8-43) 

y 

This  can  be  easily  verified  by  direct  calculation:  there  is 

2a^v  = A^<f  + y.|y-  Ay<^  . (8-44) 

The  interpretation  of  v is  as  follows.  Consider  aW  as  given  by 
(8-21).  (Of  course,  a has  here  nothing  to  do  with  Laplace's 
operator!)  It  may  also  be  expressed  in  terms  of  t by 

*w  = . . (8-45) 

In  (8-21)  , aW  has  been  considered  as  defined  on  the  telluroid  I 
It  may,  however,  also  be  regarded  as  a spatial  function,  defined 
outside  and  on  , since  t is  a function  of  the  which  can 

be  interpreted  as  curvilinear  coordinates  in  space.  If  now  aW  , 
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regarded  as  a spatial  function,  is  expressed  in  terms  of  we 

can  also  transform  (8-45)  to  these  coordinates.  This  is  best  done 
by  transforming  it  first  to  the  form 


AW  = - T . 


using  (7-29)  with  instead  of  g^  . Now 


(8-46) 


(8-47) 


is  related  to 


(8-48) 


Y = , y = /y'  . 


(8-49) 


by  (8-39).  Therefore, 


li  = y21l  1 y ll 

3y  9y  iy  ? 3y  * 


and  (8-46)  takes  the  form 


.W  . I y , 


Substituting 


T = y(j) 


(8-50) 


(8-51) 


(8-52) 


according  to  (8-40),  we  get 


AW  = 1 


1 


3y 


y*)* 
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7 y(y|y  - <(>)  • 


Using  again  {J -1^)  with  instead  of  g^^  we  obtain 


(8-53) 


1 


AW  = j yCy^f^  - <(>)  . 


(8-54) 


and  the  comparison  with  (8-42)  shows  that 


V = 


AW 


(8-55) 


This  furnishes  the  desired  physical  interpretation  of  v 
These  two  auxiliary  functions 


T_ 

y 


T 

y 


v(yi)  = 7(yi|^  - 0) 


AM 

y 


(8-56) 

(8-57) 


which,  in  the  case  of  spherical  symmetry,  satisfy  Laplace's  equation 


A (J.  = 0 , 

y 

A V = 0 , 

y 


(8-58) 

(8-59) 


will  play  a basic  role  in  the  next  section. 
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9 . The  Nonlinear  Problem 

Reformulation  of  the  Problem If  t is  a solution  of  the 
boundary-value  problem  defined  by  the  differential  equation  (7-26) 
and  the  boundary  condition  (7-31),  which,  in  view  of  (7-29),  may  be 
written  in  the  form 

0"  Sg  , (9-1) 

then  the  function 

4-  = + c.g.  , (9-2) 

with  an  arbitrary  constant  vector  c^  , is  also  a solution  of  the 
problem.  In  fact. 


ij  sgiSgj  ij 


(9-3) 


so  that  4'  satisfies  (7-26)  if  4*  does,  and 


9k 


5 


(9-4) 


so  that  the  boundary  condition  is  also  satisfied. 

It  is  easily  seen  that  the  addition  of  the  term  c^g^  to 
4"  represents  a translation  by  the  vector  c^  in  ordinary  space: 
by  (7-13)  we  get 


= Xi  + Ci 


+ c 


(9-5) 


Ill 


We  obtain  a unique  solution  by  requesting  t to  have  the  form  (7-38), 
which  places  the  x-coordinate  system  at  the  earth's  center  of  mass. 
This  is  in  complete  correspondence  with  the  usual  treatment  of 
Molodensky's  problem. 

However,  the  solution  will  not  exist  for  arbitrary  boundary 
values  V but  only  for  those  functions  T(u)  which  satisfy  n 
conditions;  from  the  discussion  of  the  linearized  problem  we  expect 
n = 3 , It  is  true  that  if  we  had  idealized  conditions,  especially 

absence  of  measuring  errors,  then  the  data  function  Y(u)  would 
satisfy  these  conditions  because  the  solution  exists  for  physical 
reasons.  In  practice,  however,  especially  because  of  measuring  and 
interpolation  errors,  we  cannot  expect  that  the  actual  V’(u)  will 
exactly  satisfy  these  conditions. 

This  suggests  a reformulation  of  the  boundary-value  problem 
in  gravity  space  along  the  lines  of  Hbrmander's  formulation;  cf.sec.6, 
especially  eq.  (6-11):  we  replace  the  boundary  condition  (9-1)  by 

9k  Ig"  ■ ^ 

The  new  boundary-val ue  problem  can  now  be  expected  to  have  a solution 
for  arbitrary  data  functions  7(0)  . The  three  constants  a^, 
are  determined  as  unknowns  and,  so  to  speak,  take  care  of  the  three 
conditions . 

T ransf ormat i on  of  the  Differential  Equation.-  The  main 
difficulty  in  the  gravity  space  approach  lies  in  the  differential 
equation  for  the  adjoint  potential  4'  . This  equation,  given  by 

(7-26)  or  (7-27),  is  unfortunately  considerably  more  complicated 
than  Laplace's  equation  for  the  original  potential  V 

The  consideration  of  the  spherical  approximation  in  the 
preceding  section  suggests,  however,  that  it  may  be  possible  to 


! 

I 

■] 


1 

\ 
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reduce,  at  least  approximately,  this  differential  equation  to 
Laplace's  equation. 

First,  in  agreement  with  (7-38),  we  split  off  the  main  part 
in  I*  by  putting 

= - 2y2g2  + ,(g_)  . (9.7) 

where 

= g^g,^  • (9-8) 

This  may  be  interpreted  by  (8-10)  as  using  a spherically  symmetric 
reference  potential 

$ = - 2y2g2  (9-9) 

in  gravity  space  (the  reader  will  find  it  best  to  consider  all 
transformations  to  follow  as  transformations  in  gravity  space  and 
to  forget,  for  the  time  being,  about  ordinary  space).  In  contrast 
to  the  linear  treatment  in  the  preceding  section  we  shall  not  intro- 
duce any  approximations,  so  that  the  transformed  differential 
equations  will  be  as  rigorous  as  the  original  one,  eq.  (7-26). 

The  reference  potential  (9-9),  which  is  spherically  symmetric 
in  gravity  space,  is  the  adjoint  potential  of  a potential  V that 
is  spherically  symmetric  in  ordinary  space.  In  fact,  by  (7-13)  and 
(7-15), 


113 


I 2 ~ 2 

r = = M g 


(9-11) 


1 1 


1 1 


V = 


9k^k 


2„2 

p g 


y = 


£ 

r 


2„2  . , 2„2 

M g + 2m  g 


(9-12) 


Eq.  (8-31)  suggests  the  substitution 


= g 


9i 


(9-13) 


(Now,  however,  the  are  to  be  considered  as  curvilinear  coordi 

nates  in  gravity  space,  having  no  direct  relation  with  cartesian 
coordinates  in  ordinary  space.)  This  transforms  the  reference 
potential  into 


y = - 2m  y 


(9-14) 


eliminating  the  singularity  g at  the  origin 
We  now  introduce  the  new  function 


T_ 

y 


(9-15) 


so  that 


yi* 


(9-16) 


If  we  neglected  all  squares  and  higher  powers  of  t , we  should 
have  the  linear  spherical  approximation  discussed  in  the  preceding 
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section  since,  apart  from  a scale  factor,  (9-12)  is  identical  to 
(8-27).  This  shows  that  t , as  a function  of  y , must  satisfy' 
a differential  equation  of  form 


= 0((^  ) ; 


(9-17) 


there  can  be  no  term  0(4,)  on  the  right-hand  side  since  a^(()  = 0 
as  a linear  approximation,  by  (8-58). 

In  fact,  Sansd  has  calculated  the  exact  differential  equa- 
tion which  (()  must  satisfy.  This  is  done  by  substituting 


f = - 2y^y  + y.^ 


(9-18) 


into  (7-27)  and  performing  some  lengthy  but  straightforward  trans- 
formations. The  result  is  ( Sanso , 1977a  , p .69 ) : 


A * = M (4)  ,<)> ) , 


where  8^(4), 4))  is  a quadratic  operator  given  by 


(9-19) 


B^(0>.*)  = ^(<0  - yi.')A^4.  + y2  [(TrL)2  - Tr(L2)]  (9-20) 

(it  must  be  quadratic  since  the  original  equation  (7-26)  is).  The 
matrix  L has  elements 


L . , = (6.^ 

1 j ' ik 


T — 


(9-21) 


where  6 


denotes  the  elements  of  the  unit  matrix  and 


I 


f f 
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(|) ' is 


d,  = _il_2 

j-d  9y^9yj 

defined  by 


(9-22) 


and  A^<p  expresses  the  Laplace  operator  in  the  "cartesian"  form 

A ^ + ill  + . (9-23) 

^ 3yj 

needless  to  say,  y^  are  not  rigorously  to  be  interpreted  as 
cartesian  coordinates  in  ordinary  space. 

The  boundary  operator  (9-1) 


, - f = g 21  . ^ 

'k  ag,  y ag  ^ 


is  transformed  as  follows.  Using  (9-9)  we  find 

„ ay  - _ 2„2 

9 “ 9 • 


SO  that 


3'*'  - 2„2  . „ dT 

ag  ^ ^ ^ ag 


Since 


(9-24) 


(9-25) 


(9-26)  I 

j 

( 

i 


(9-27) 
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by  ( 9- 13 ) , we  have 


ll  = 

3g  3y  "ag  7^  3y 


In  view  of  these  relations  we  get 


1 IL 

k 39, 


't'  = y^y 


iy—  ' 

3y 


(9-28) 


(9-29) 


On  substituting  (9-16)  and  taking  (9-6)  into  account  we  find  as 
boundary  condition  for  4)  : 


y 


9 4) 

9y 


v(u)  + 


a.y  1 

i-^  i I 


on 


(9-30) 


where 


1 


(9-31) 


is  a function  of  the  data,  and 

y(u)  = yoS^  , y^(u)  = y^oS^  (9-32) 

denote  the  values  of  y and  y.  calculated  for  that  point  of  the 
surface  S which  has  the  parameter  u 

Since  the  direction  of  3/3y  , as  well  as  the  direction  of 

3/39  » is  the  direction  of  the  radius  vector  in  gravity  space, 

we  still  have  an  oblique  derivative  problem  as  in  the  ori9inal 
formulation  given  at  the  end  of  sec. 7;  the  problem  is,  however. 
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simplified  because  we  now  have  a "quas i 1 i near " differential 
equation  (9-19),  which  has  a form  suitable  for  an  iterative  solution. 

It  is,  nevertheless,  appropriate  to  transform  the  problem 
still  further  by  introducing  a new  potential  v by 

' ' 

This  substitution  has  been  motivated  in  the  preceding  section;  cf. 
(8-57).  As  a linear  approximation,  v(y^)  is  harmonic  and  is, 
furthermore,  related  to  the  potential  anomaly  AW 
In  fact,  we  have  even  rigorously 

v(g.)  = + yv  , (9-34) 

so  that  yv  represents  the  perturbation  in  the  potential  V , if 
expressed  in  gravimetric  coordinates  g^  , in  the  same  way  as  we 
had 


v(gj  = - + y^  , (9-35) 

T = y(t)  representing  the  perturbation  in  the  adjoint  potential  4' 
It  is  easy  to  verify  (9-34)  by  substituting  (9-35)  into  (7-30). 

By  means  of  (9-33),  eq.  (9-19)  is  finally  transformed  into 
a differential  equation  for  v 

AyV  = y ^B2(v,v) 


(9-36) 
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where  the  quadratic  operator  82  is  given  by 

y 

B_ (v,v)  = - VA  V - V ' JA  vdy  + 

I y ' y 

O 


+ 2 [(TrN)^  - Tr(N^)] 

+ 4y 

TrM'TrN  - Tr(MN) 

(9-37) 

M and 

N 

are  3x3  matrices  with  elements 

M 

- 

(9-38) 

N 

i j 

y 0 

(9-39) 

where 

6 . 
3-: 

denotes  the  elements  of  the  unit  matrix  and 

V . , 

3-: 

2 

3 V 

» 

(9-40) 

V ' is 

defined  as 

V ' 

_ 3V 

ay 

(9-41) 

and  A^  again  denotes  the  "cartesian  form"  of  the  Laplace  operator. 

The  solution  of  (9-36)  has  to  satisfy  the  boundary  condition 

voS  = V + ay 


(9-42) 
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V = 7(u)  and  y = y(u)  being  given  by  (9-31)  and  (9-32).  This  is 
simply  the  boundary  condition  for  a Dirichlet  problem. 

By  means  of  the  substitution  (9-33)  it  has  thus  been 
possible  to  transform  Sansd's  problem  into  a Dirichlet  problem 
for  the  nonlinear  equation  (9-36).  The  price  to  be  paid  is  that  this 
equation  is  a nonlinear  i ntegro-di ff erential  equation,  as  (9-37) 
shows.  However,  since  the  principal  part  of  (9-36)  is  simply  Laplace's 
equation,  the  quadratic  right-hand'  side  being  relatively  small,  our 
equation  is  still  relatively  manageable  (it  is  hardly  necessary  to 
remind  the  reader  that  (9-36)  is  as  rigorous  as  the  original  equation 
(7-26),  no  neglections  are  involved). 

This  reduction  to  a Dirichlet  problem  is  similar  to  methods 
used  in  the  linear  Molodensky  problem,  cf.  (Brovar , 1964) , (Krarup, 
1973:  the  "Prague  method")  and  the  reduction  the  Brillouin  sphere 
in  the  present  sec. 4.  The  enormous  advantage  of  the  gravity  space 
approach  is  that  the  boundary  condition  (9-1)  is  linear  even  for  the 
nonlinear  problem,  so  that  methods  can  be  used  that  are  applicable 
to  the  Molodensky  problem  only  in  its  linearized  form. 

A necessary  condition  for  the  existence  of  the  solution  is 


(9-43) 


In  fact,  the  differentiation  of  (9-33)  gives 


2 av  ^ 


- V X 


(9-44) 


If  the  solution  iji  is  to  be  regular  with  finite  second  derivatives 
at  the  origin  y = /y^y ^ = 0 , then  (9-44)  must  tend  to  zero  as 
y^  0 . The  condition  (9-43)  is  to  be  provided  for  by*suitably 
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disposing  of  the  free  constants  a^  in  the  boundary  condi 

tion  (9-30) . 

If  a solution  v satisfying  (9-43)  has  been  found,  then 
41  is  obtained  by 


4.  = - 2v(0)  + 


v(0)]  y~^dy 


(9-45) 


it  is  easy  to  verify  by  direct  substitution  that  this  solution 
satisfies  (9-33).  Then  the  adjoint  potential  4'  is  given  by  (9-35), 
and  finally  the  earth's  surface  is  obtained  by  (7-39).  A check  is 
provided  by  (9-34) . 

Compared  with  the  Nash-Hbrmander  approach  to  the  nonlinear 
Mol  Odens ky  pro bl em,  it  is  relatively  simple  and  straightforward  to 
solve  (9-36)  by  a Newton  iteration  scheme,  considering  p'^'^^as  a 
small  parameter.  In  this  way,  Sansd  has  obtained  first  results  on 
existence  and  uniqueness  of  the  solution.  He  has  proved  that  a 
unique  solution  exists  provided 

< 6 , (9-46) 

where  the  constant  6 is  sufficiently  small.  The  norm  is  a Holder 
norm  very  similar  to  the  norm  used  in  sec. 6. 

The  condition  (9-46)  is  directly  comparable  to  (6-25);  the 
constants  6 will  be  different.  It  means  that  V as  given  by  (9-31) 
should  be  small,  as  well  as  the  first  and  second  derivatives,  inclu- 
ding a Holder  condition  on  the  second  derivatives.  Since  v denotes 
the  deviations  from  a spherical  symmetrical  solution,  this  condition 
is  again  very  restrictive. 


r 
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10.  Conclusions 


During  the  recent  years,  the  problem  of  existence  and 
uniqueness  of  the  solution  for  Molodensky's  problem  has  for  the 
first  time  been  treated  with  adequate  mathematical  rigor.  Certainly, 
existence  and  uniqueness  have  been  proved  only  under  very  restrictive 
conditions  on  smoothness  and  smallness  of  the  deviations  from  a 
"normal"  solution,  conditions  which  are  hardly  met  in  the  actual 
geodetic  situation.  However,  these  results  have  been  obtained 
rigorously. 

The  treatment  by  Hbrmander  uses  a very  advanced  inverse 
function  theorem  and  is  mathematically  extremely  complicated;  it 
applies  to  a rotating  earth.  The  mathematical  complexity  is  mainly 
due  to  the  fact  that  Molodensky's  problem  is  a free  boundary-value 
problem,  the  boundary  surface  being  unknown. 

The  gravity  space  approach  due  to  Sansd  transforms  the  free 
boundary  problem  into  a fixed  one,  although  for  a nonlinear  partial 
differential  equation.  It  nevertheless  reduces  essentially  the 
mathematical  complexity.  The  limitation  of  the  gravity  space  approach 
is  the  restriction  to  a nonrotating  earth;  practically  this  amounts 
to  the  use  of  gravitation  instead  of  gravity  by  reducing  for  the 
effect  of  centrifugal  force. 

From  this  point  of  view,  an  extension  of  Sansd's  approach 
to  a rotating  earth  by  an  iterative  procedure  (using  the  fact  that 
u)  is  small)  appears  less  urgent;  at  any  rate  it  seems  to  be  not 
quite  easy  to  prove  convergence  of  such  an  iteration. 

The  results  obtained  so  far  by  Hbrmander  for  co  ^ 0 and 
by  Sansd  for  w = 0 are  comparable;  the  conditions  are  similarly 
restrictive.  It  would,  of  course,  be  desirable  to  obtain  stronger 
results,  for  norms  |1  or  even  better  |1  1|^  with  0<e<l 
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For  this  purpose  the  gravity  space  approach  appears  to  be  more 
promising  since  it  is  so  much  easier. 

The  impact  of  the  gravity  space  approach  to  the  theory  of 
Molodensky's  problem  appears  to  be  enormous;  it  may  well  be  compa- 
rable to  the  impact  of  Hamiltonian  methods  to  Newtonian  classical 
mechanics  (both  apply  a Legendre  transformation!). 

From  a practical  point  of  view  it  is  important  to  note  that 
the  linear  approximation  (linear  in  the  anomalous  potential  T ), 
which  is  probably  sufficient  for  almost  all  present  purposes,  is  the 
same  in  the  usual  approach  and  in  gravity  space.  Therefore,  the 
number  of  the  usual  methods  for  practically  solving  Molodensky's 
problem,  such  as  Molodensky's  series  and  related  solutions,  is  not 
augmented  by  the  new  developments. 
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